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Eastman  Kedah  Gempany  recently  completed  contract  with  the  United  States 

Air  Force  Air  Weather  Service  to  conduct  research  into  the  transfer  of 
atmospheric  energy  in  the  visible  and  near- infrared  portions  of  the  elec- 
tromagnetic spectrum.  Its  objective  was  to  produce  a spectral  model  of 
path  transmittance,  radiance,  and  ground- level  irradiance  that  could  be  re- 
lated to  meteorological  observations  through  the  use  of  simultaneous,  in- 
situ  data  collections.  For  this  purpose  a mathematical  model  was  developed. 
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infrared  spectral  absorption  and  scattering  phenomena  that  result  from  atmo- 
spheric constituents,  from  surface  synoptic  weather  observations. 
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radiometric  and  meterological  data  to  a series  of  equations  which  became  the 
SCAT3  model.  Examples  are  also  shown  illustrating  the  final  model's  ability 
to  match  observed  spectral  structure  based  upon  coefficients  derived  from 
surface  meteorology. 
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SECTION  1 
INTRODUCTION 


The  approach  used  to  develop  the  equations  of  the  SCAT3  atmospheric 
model  was  to  start  with  a single  linear  mathematical  description  of  the 
physical  processes.  The  results  of  a detailed  examination  of  changes 
in  the  spectral  data  were  then  used  to  introduce  non-linear  effects  into 
coefficients  of  the  various  linear  terms.  The  net  effect  of  this  approach 
was  to  create  parametric  spectral  vectors  that  could  be  fitted  by  least- 
squares  to  appropriately  normalized,  measured  spectral  data.  The  coeffi- 
cients of  the  vectors  could  then  be  regressed  against  other  measurable 
parameters  such  as  solar  altitude,  moisture  scale  height,  and  terrain  re- 
flectance. In  this  way,  a model  could  be  assembled  that  closely  approx- 
imated the  physics  of  the  atmosphere,  thus  avoiding  the  equation  extrapo- 
lation problems  that  are  often  encountered  in  a characteristic  vector  ap- 
proach to  data  modeling. 

This  approach  was  used  for  all  radiometric  data,  which  included  fl)  hori- 
zontal daylight  and  skylight  irradiance,  (2)  vertical  front  daylight  and 
vertical  back  skylight  irradiance,  and  (3)  atmospheric  path  radiance  and 
transmittance.*  In  all  the  cases  except  transmittance,  the  measured  spec- 
tral data  was  normalized  by  the  extra-terrestrial  solar  irradiance.  The 
mathematical  quantities  representing  direct  beam  attenuation,  primary 
atmospheric  scattering,  and  surface  reflectance  were  then  entered  into  the 
equation.  Regression  coefficients  were  computed  by  fitting,  by  least-squares, 
the  spectral  structure  to  each  normalized  spectral  sample.  The  coefficients 
were  then  analyzed  as  functions  of  such  measurable  parameters  as  solar  alti- 
tude, moisture  scale  height,  and  terrain  reflectance. 

Data  was  blocked  in  small  increments  of  solar  altitude.  This  was  possible 
because  of  the  large  volume  of  spectral  data  the  mobile  radiometric  labora- 
tory was  able  to  collect.  The  size  of  solar  altitude  cells  into  which  the 


* See  Appendix  A for  a description  of  the  software  developed  for  this 
purpose. 
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data  was  "blocked  (actually  interpolated  to  a common  set  of  solar  altitudes) 
varied  with  solar  altitude;  that  is,  between  -3°  and  10°  the  increment 
was  1°,  between  10°  and  50°  it  was  2°,  and  between  50°  and  90°  it  was  5°. 

The  increment  chosen  was  quasi- logarithmic  and  corresponuea  roughly  to 
equal  increments  of  change  in  the  measured  variable.  The  effects  of  the 
remaining  parameters,  moisture  scale  height  and  terrain  reflectance,  were 
determined  sequentially;  this  was  done  by  examining  least-squares  fits  to 
the  residual  spectral  structure  that  were  obtained  when  predetermined  equa- 
tional  shapes  were  regressed  against  the  actual  data.  These  equations  were 
based  on  approximations  of  the  integrated  effects  of  the  atmosphere  without 
recourse  to  a multilayered  atmosphere;  in  this  way,  the  formulation  of  terms 
was  an  approximation  of  the  physical  equations.  Several  alternative  equational 
forms  were  examined,  but  the  ones  used  in  this  effort  consistently  resulted 
in  the  best  fit  to  the  spectral  structure. 


The  regression  coefficients  corresponding  to  scattering  and  surround  reflec- 
tance effects  were  unique  for  each  measured  spectral  sample.  (For  an  example, 
see  the  equation  in  para  4.2.)  However,  by  studying  first  the  samples  for 
which  the  surround  term  could  reasonably  be  expected  to  be  zero,  it  was  possible 
to  determine  the  scattering  relationship  to  moisture  scale  height.  Once  the 
scattering  function  was  determined,  the  surround  relationship  to  terrain  re- 
flectance could  be  found.  Because  these  terms  were  not  totally  independent 
of  one  another  in  practice,  some  complication  was  introduced  into  the  analysis. 


Similar  equations  were  built  up  for  all  four  types  of  spectral  irradiance  and 
spectral  sky  radiance.  The  resulting  formulation  of  each  can  be  found  in 
appropriate  sections  of  this  volume.* 

The  ability  of  the  approach  used  here  to  separate-out  effects  necessitated 
having  precise  estimates  of  spectral  extra-terrestrial  irradiance,  spectral 
transmittance,  and  spectral  terrain  reflectance.  The  degree  of  accuracy  re- 
quired was  determined  by  each  term's  influence  on  calculation  of  the  spectral 


*See  Appendix  B for  a summary  of  the  variable  names  used  in  the  equations  and 
in  the  SCAT3  model. 
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radiometric  quantity.  The  extra-terrestrial  irradiance  selected  for 
normalization  of  the  spectral  samples  was  the  so-called  Johnson  Solar 
Constant1  corrected  for  the  earth's  orbital  eccentricity.  The  solar  con- 
stants discussed  by  Neckel2  and  Thekaekara3  were  also  examined,  but  best 
agreement  with  measured  data  was  generally  obtained  with  the  Johnson  con- 
stant. 

Estimates  of  the  spectral  transmittance  due  to  scattering  and  absorption, 
t (A)  and  t (A) , were  obtained  simultaneously  from  direct  measurement  of  the 
solar  disk  radiance  through  the  atmosphere.  The  procedure  followed  in  ana- 
lyzing and  decomposing  spectral  transmittance  into  its  component  parts  is 
described  in  Section  3 of  this  volume. 

Estimates  of  spectral  terrain  reflectances  for  each  location  date  were  de- 
rived from  color  overhead  photography,  or  in  some  instances,  from  geo-botanical 
maps.  By  means  of  a spectral  Eigenvector  technique  and  knowledge  of  the  type 
and  relative  amounts  of  groundcover,  it  was  possible  to  expand  the  three  inte- 
grated reflectances  derived  from  analysis  of  the  color  film  dye  layer  into  a 
complete  spectral  reflectance  estimate  for  the  general  terrain  reflectance 
corresponding  to  a 50-  to  100-mile  area  surrounding  the  site  of  operation. 

These  estimates  are  discussed  in  Section  S of  this  volume. 

A matrix  of  actual  measured  spectral  samples  was  created  that  represented 
typical  and  extreme  values  of  model  parameters  and  were  descriptive  of  the 
standard  error  of  spectral  reconstruction.  These  reconstructed  samples  are 
shown  in  Appendix  C of  this  volume.  Because  of  the  importance  of  analysis  of 
spectral  transmittance  to  the  results  of  the  model's  fit  to  data,  this  subject 
is  discussed  below  before  the  topics  spectral  radiance  and  irradiance. 


1 F.S.  Johnson,  "The  Solar  Constant,"  J.  Meteorology,  Vol  II,  No.  6, 

December  1954. 

2 D.  Neckel,  "Intensity  of  the  Center  of  the  Solar  Disk  in  the  Spectral 
Region  .33  to  1.25  run  Measured  from  High  Mountain  Station,"  Proceedings 
of  Symposium  on  Solar  Radiation,  November  1973,  Smithsonian  Institute. 

3 Dr.  M.  P.  Thekaekara,  "Solar  Electromagnetic  Radiation,"  NASA  Space  Vehicle 
Design  Criteria  (Environment) , NASA  Sp.  83005, Rev  May  71,  Goddard  Space 
Flight  Center. 
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SECTION  2 


ANALYSIS  AND  MODELING  OF  SPECTRAL  TRANSMITTANCE 

2.1  APPROACH 

Atmospheric  beam-path  spectral  transmittance  was  used  as  the  basic  scal- 
ing parameter  of  all  the  spectral  radiometric  quantities  computed  by 
SCATS.  As  such,  it  was  the  single  most  important  measurement  made  by  the 
mobile  laboratory  and  the  most  critical  in  the  regression  analyses  that 
are  discussed  in  the  following  sections.  Spectral  transmittance  also  pro- 
vided a link  between  the  radiometric  quantities  and  the  detailed  altitude 
profile  data  on  the  atmosphere  provided  by  the  Rawindsonde  observations 
made  at  the  site  of  operations. 

The  approach  taken  in  the  regression  analysis  was  to  decompose  the  loga- 
rithm (base  e)  of  the  spectral  transmittance  into  its  individual  spectra 
in  the  visible  and  near-infrared.  The  assumption  was  made  that  molecular 
scattering  transmittance  could  be  approximated  by  the  Rayleigh  equations 
scaled  by  atmospheric  pressure.  Aerosol  scattering  transmittance  was 
estimated  as  the  residual  transmittance  between  the  measurement  and  the 
Rayleigh  equations,  eliminating  the  absorption  transmittance  wavelengths. 
Absorption  transmittance  was  then  estimated  by  fitting  the  linear  charac- 
teristic vectors  to  the  measured  data  after  the  scattering  effects  were 
subtracted. 

Meteorological  data  were  acquired  as  a part  of  this  project  for  the  pur- 
pose of  examining  the  correlations  with  atmospheric  transmittance  that  had 
been  discussed  by  Harrel  and  Bullrich4and  other  researchers  in  the  field  of 
atmospheric  science.  In  the  study,  measurements  were  made  of  temperature, 
pressure,  relative  humidity,  wind  speed,  and  wind  direction  using  an  on- 
board weather  system.  In  addition,  Rawindsonde  (RAOB)  observations  that 


4 G.  Harrel  and  K.  Bullrich,  "On  the  Interpretation  of  Atmospheric  Turbidity 
Measurements,"  Journal  of  the  Atmospheric  Sciences,  Vol  S3,  pp  794-797 
May  1975. 
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were  made  three  times  a day  by  the  National  Weather  Service  were  used 
as  the  source  of  altitude  profile  data  on  temperature,  pressure,  and  mois- 
ture. The  NWS  surface  data  were  also  collected,  where  available,  and  they 
were  given  preference  over  on-board  weather  measurements.  Correlations 
were  found  to  exist  between  meteorological  data  and  atmospheric  transmit- 
tance both  in  the  wavelength  bands  corresponding  to  water  absorption  bands 
and  in  the  aerosol  scattering  bands. 

The  schematic  diagram  in  Figure  2-1  shows  the  physical  relationship  be- 
tween the  various  data  collected  by  the  mobile  laboratory. 

2.2  ESTIMATION  OF  VERTICAL  TRANSMITTANCE  FROM  MEASUREMENTS 

Direct  estimation  of  the  solar  disk  radiance  inside  the  atmosphere  was 
possible  because  of  the  on-board,  absolute  calibration  of  the  transmis- 
someter  (see  Vol  1,  para.  3,1  of  this  report).  The  calibration  was  carried 
out  in  a manner  that  obtained  data  directly  from  the  instrument  current 
output  in  units  of  path  transmittance.  After  careful  study,  it  was  deter- 
mined that  on  an  absolute  scale,  the  instrument  demonstrated  unacceptable 
non-linearities  in  the  blue-UV  spectrum  that  resulted  in  underestimation 
of  the  path  transmittance  in  that  region.  The  balance  of  the  spectrum, 
from  0.5  pm  to  1.15  pm,  appeared  to  be  in  agreement  with  theory  and  with 
alternative  methods  of  transmittance  estimation. 

The  alternative  method  that  was  examined  and  ultimately  adopted  for  the 
analysis  of  spectral  radiometric  data  discussed  in  Section  1 of  this  volume 
was  a modification  of  the  Langley-plot  method  of  solar  disk  radiance  deter- 
mination. This  method,  which  was  originally  intended  for  the  estimation  of 
extra-trrestrial  solar  irradiance,  involves  plotting  the  logarithm  of  disk 
radiance  versus  the  relative  air  mass  of  observation  as  estimated  using  the 
secant  of  the  observation  zenith  angle  or  air  mass  tables  developed  by 
Bemporad5  at  low  solar  altitudes.  The  intercept  of  a straight  line  fitted 


5 A.  Bemporad,  "Search  for  a New  Empirical  Formula  for  the  Representation 
of  the  Variation  of  the  Intensity  of  Solar  Radiation  with  Zenith  Angle," 
Meteorologische  Zeitschrift,  Vol  24,  July  1907. 
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to  the  radiance/air  mass  plot  at  zero  air  muss  was  the  extrapolated 
estimate  of  the  extra-terrestrial  radiance.  Normalization  of  the  radiance 
scale  by  the  solar  constant  resulted  in  a path  transmittance  versus  air 
mass  that  was  consistent  with  the  calibrated  output  of  the  transmissometer 
instrument.  Ideally,  the  linear  fit  to  the  data  would  result  in  an  inter- 
cept of  0.0  (log  1.);  however,  the  extrapolated  estimation  was  generally 
less  than  zero,  particularly  below  a wavelength  of  500  nm. 

Examination  of  the  equation  of  the  slope  of  log^  of  the  transmittance;  that 
is,  log  T(X)  versus  the  relative  air  mass  m,  if  it  is  assumed  linear,  yields 
an  instrument  calibration  estimate  independent  of  the  vertical  atmospheric 
transmittance,  as  follows: 

T(X)  = log  Tv(X)  Equation  2-1 

Using  a simple  linear  regression  technique,  slope  estimates  of  the  vertical 
transmittance  were  derived  for  every  wavelength  and  every  site  both  before 
or  after  true  noon.  A typical  example  of  such  a plot  for  three  wavelengths 
at  one  site  is  shown  in  Figure  2-2  for  Vandenberg,  California,  on  May  12, 

1976.  Not  all  day-sites  exhibited  a perfectly  linear  relationship,  but  by- 
and- large,  all  scattering  wavelength  data  could  be  approximated  by  a linear 
function  within  the  standard  error  of  the  measurement  (estimated  to  be  less 
than  5 percent,  see  Vol  1,  para.  3.4).  Absorption  transmittance  measure- 
ments corresponding  to  the  water  and  molecular  oxygen  absorption  bands  also 
exhibited  non-linear  variation  with  air  mass  because  these  absorption  spectra 
are  not  strictly  linear  with  total  absorber  concentration,  as  discussed  by 
McClatchey.6  The  analysis  proceeded  linearly,  however,  with  a correction 
applied  in  a later  stage  of  analysis  for  the  non-linearity  of  the  observed 
data. 


6 R.  A.  McClatchey  et.  al.,  Optical  Properties  of  the  Atmosphere  (Third 
Edition) , Air  Force  Cambridge  Research  Laboratory  Project  7670. 


2.3  ANALYSIS  OF  VERTICAL  SPECTRAL  TRANSMITTANCE  ESTIMATES 


By  scaling  samples  linearly  for  air  mass  and  correcting  for  the  absorption 
nonlinearity,  any  individual  spectral  transmittance  sample  could  be  re- 
constructed. These  reconstructed  samples  for  each  day-site  were  the  basic 
spectral  transmittance  data  used  in  all  subsequent  radiometric  analysis. 

An  example  of  the  measured  spectral  transmittance  showing  the  spectral 
structure  and  using  10-nanometer  instrument  sampling  is  shown  in  Figure 
2-3;  this  measured  spectral  transmittance  was  estimated  by  the  modified 
Langliy-plot  method. 


The  equivalent  spectral  transmittance  for  each  day-site  was  then  logged 
(to  the  base  e)  and  decomposed  into  Rayleigh  and  Mie  scattering  terms  and 
into  molecular  absorption  spectra.  This  was  done  in  a series  of  steps  begin- 
ning with  the  removal  of  an  assumed,  fixed  amount  of  total  path  transmittance 
due  to  Rayleigh  scattering  that  was  scaled  for  atmospheric  surface  pressure 
in  millibars  as  described  by  Robinson.7  If  the  absorption  bands  are  tempor- 
arily ignored,  the  net  spectra' is  that  corresponding  to  total  extinction  due 
to  large  particle  scattering.  It  was  hypothesized  that  this  scattering  re- 
sulted from  a combination  of  dust  particles  and  water  droplets  that  com- 
prise the  so-called  atmospheric  aerosols.  It  was  further  hypothesized  that 
this  would  scale  with  the  water  vapor  content  of  the  atmosphere.  A non- 
linear function  of  the  form  suggested  by  Robinson7  and  shown  below  was  fitted 
to  the  net  spectra  using  regression  techniques. 


where  Kft  is  the  aerosol  scattering  coefficient  and  c is  the  aerosol  scat- 
tering wavelength  exponent.  Figure  2-4  illustrates  the  fit  of  the  function 


7 N.  Robinson,  Solar  Radiation,  Elsevier  Publishing  Co.,  New  York, 
1966,  p 113. 
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Figure  2-3.  Example  of  Measured  Spectral  Vertical  Atmospheric  Log 
Transmittance  (Vandenberg  , Ca.;  May  12,  1976  P.M.) 
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Figure  2-4.  Typical  Regression  Model  Fitted  to  Net  Aerosol 
Spectral  Vertical  Atmospheric  Log  Transmittance 
(Vandenberg,  Ca.;  May  12,  1976  P.M.) 
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to  the  net  log  transmittance  data.  Recombining  the  scattering  phenom- 
ena equations  and  subtracting  them  from  the  total  measured  log  transmit- 
tance resulted  in  the  log  transmittance  spectra  for  absorption  alone. 

The  equation  for  total  absorption  extinction  was  then: 


l0*W» 


logT  (A)  - K /A4  - K /AC 
in  r a 


Equation  2-2 


As  was  expected,  except  for  the  spectral  regions  affected  by  absorp- 
tion, the  net  spectra  is  zero  (log  E 1.0)  within  the  standard  error  of  the 
equation's  fit  to  the  measured  data.  Arbitrary  scale  height  values  were 
then  found  by  linear  regression  techniques  for  molecular  ozone,  oxygen, 
and  water,  all  of  which  display  sizable  absorption  spectra  in  this  wave- 
length region.  These  bands  are  clearly  visible  in  Figure  2-5,  which  is 
an  example  of  the  model  equation  used  to  approximate  all  structure  of 
each  day-site  spectral  transmittance.  The  following  equation  combines 
these  terms.  Figure  2-6  shows  this  equation  fitted  to  the  total  measured 
log  transmittance  for  Vandenberg,  Ca. ; May  12,  1976  p.m. 


•Exp  - Kr/A4  - K&/ 


K V (A)  - K V (A)  - K V 
o o ' y oo 
2 2 3 


,w] 


Equation  2-3 


Where:  K ,K  , and  K are  the  relative  scale  weight  values 
w °2  °3 

* corresponding  to  water,  oxygen,  and 
ozone  absorption,  and 

V (A),  V (A),  and  V (A)  are  the  relative  spectral  shape  vectors 
w o o, 

for  water,  oxygen,  and  ozone  absorption 
spectra. 

For  the  sample  date  used  in  Figures  2-2  through  2-6,  the  values  of  the 
constants  of  each  term  were: 


Ka  • 0.0712 

c • 0.685 


0.100 


k ■ 0.0157 

°3 


K • 0.107 

°2 

K • 0.00884 
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In  subsequent  analyses,  the  relative  scale  heights  were  related  to  the 
physically  measurable  parameters  discussed  below  in  this  section.  However, 
for  purposes  of  analysis  of  spectral  sky  radiance  and  spectral  irradi- 
ance  data,  the  exact  coefficients  derived  from  this  reduction  process  were 
used.  These  coefficients  are  summarized  in  Table  2-1  for  all  day  sites 
for  which  the  transmittance  analysis  could  be  accomplished. 

2.4  MODELING  OF  VERTICAL  SPECTRAL  TRANSMITTANCE 

2.4.1  Scaling  for  Molecular  Scattering 

The  molecular  scattering  beam  alternation  can  be  adequately  described  by 
Rayleigh  scattering  theory.  The  coefficient  of  the  Rayleigh  form  has  been 
shown  to  be  a function  of  the  atmospheric  density  ana  thus  proportional  to 


-0.008925 

1013.25 


x AP 


Equation  2-4 


For  the  earth-to-space  problem,  AP  becomes  the  uncorrected  surface  pres- 
sure in  millibars.  For  observations  from  an  arbitrary  altitude,  AP  is  the 
pressure  difference  between  the  observer's  altitude  and  the  terrain  elevation, 
and  it  can  be  found  by  integrating  the  atmospheric  density  versus  the  altitude 
function.  This  density  function  has  been  modeled  in  SCAT3  as  a closed-form, 
continuous  function  of  altitude  referenced  to  sea  level.  Each  ROAB  density 
profile  was  fit  with  (1)  the  function  developed  by  Weibull8  and  (2)  the  co- 
efficients studied  as  functions  of  surface  meteorological  measurements  that 
were  made  concurrently.  The  density  function  and  its  integral  between  any  two 
elevations  are  expressed  in  the  following  two  equations: 


D(x)  = Pc 


Equation  2-5 


8 W.  Weibull,  "A  Statistical  Distribution  Function  of  Wide  Applicability," 
J.  of  Applied  Mechanics,  Vol  18,  Mar  1951,  pp  293-297. 
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TABLE  2-1 

TRANSMITTANCE  COEFFICIENTS  SUMMARIZED  BY  DAY-SITE 


r 


Dale 

Sill 

Rayleigh 

Coif  ft  dm  t 

Aerosol 

Coefficient 

“a 

Aerosol 

Exponent 

C 

Niter 
Absorption 
Coif'  'ml 

*W 

Osone 

Absorption 

Coiffidmt 

*03 

Oxygen 

Absorption 

Coefficient 

*02 

«-l w 

Athena,  Cl. 

O.OQ872 

0.176 

1.384 

a 097 

0.094 

8-27-74  pa 

Apelachioola,  Fli. 

0.00899 

0.199 

1.486 

0.173 

0.163 

0.0034 

8-28-75  s- 

Apalachicola,  Fli. 

0.00899 

0.347 

1.068 

0.189 

O.I64 

0.0108 

11-9-75  pa 

Vint  Kill,  Vi. 

0.00890 

0.0648 

1.244 

0.114 

0.041 

_ 

11-11-74  pi 

Vint  Hill,  Vi. 

0. 00892 

0.0936 

0.241 

0.078 

0.007 

_ 

10-21-74  pi 

Denver,  Col. 

0.00736 

O.O484 

0.744 

0.037 

0.013 

0.034 

2-21-76  pi 

St.  Cloud,  Minn. 

0.00848 

0.0606 

0.748 

0.037 

0.044 

0.029 

2-22-76  f» 

St.  Cloud,  Minn. 

0.00870 

0.0436 

0.676 

0.043 

0.047 

0.033 

2-24-76  aa 

St.  Cloud,  Minn. 

0.00848 

0. 0690 

0.801 

0.040 

0.034 

0.069 

2-24-76  |a 

St.  Cloud,  Minn. 

0.00846 

0.0473 

0.899 

0.043 

0.043 

0.029 

4-4-76  m 

Clisgou,  Mont. 

0.00826 

0.0483 

0.989 

0.020 

0.077 

0.033 

3-5-76  as 

Ollsgou,  Mont. 

0.00824 

0.0406 

0.742 

0.012 

0.061 

0.046 

4-26-76  as 

Alaaeda,  Cil. 

00894 

0.0410 

0.414 

0.040 

0.034 

0.142 

4-26-76  pa 

AlMWdl,  Cll. 

0.00894 

0.0740 

0.918 

0.023 

0.007 

0.121 

4-12-76  a 

Vendenbirg,  Cll. 

0.00884 

0.0429 

1.130 

0.112 

0.0072 

0.130 

4-12-76  pa 

Vandmbirg,  Cal. 

0.00884 

0.0712 

0.684 

0.100 

0.016 

0.107 

6-7-76  pa 

Tucion,  Aril, 

0.00810 

0.0390 

0.190 

0.042 

0.029 

0.086 

6-6-76  pa 

Tucion,  Aril. 

0.00811 

0.0340 

0.648 

0.064 

0.023 

0.173 

6-11-70  mi 

Tucion,  Aril. 

0.00612 

0.0490 

0.184 

0.070 

0.022 

0.089 

6-17-76  m 

Tucion,  Arli. 

0.00811 

0.0381 

0.792 

0.080 

0.023 

0.086 

6-28-76  mi 

Wlnclov,  Aril. 

0.00742 

0.0639 

1.192 

0.083 

0.043 

0.157 

7-19-76  M 

Llttll  Rock,  Ark. 

0.00889 

0.062R 

1.747 

0.183 

0.014 

0.134 

7-21-76  aa 

Ldttll  Rock,  Ark. 

0.00884 

0.294 

0.342 

0.112 

0.018 

0.100 

6-4-76  ai 

AUimi,  Oi. 

0.00876 

0.1311 

1.263 

0.119 

0.078 

0.083 

6-6-76  as 

Athmi,  Cl, 

0.00872 

0.203 

0.997 

0.108 

0.166 

0.124 

6-11-76  mi 

Athmi,  Oi. 

0.00877 

0.143 

1.247 

0.091 

0.081 

0.094 

6-12-76  mi 

Athens,  Oi. 

0.00877 

0.184 

1.304 

0.143 

0.072 

0.091( 

6-13-76  aa 

Athmi,  Ci. 

0.00871 

0.209 

1.124 

0.148 

0.122 

0.110 

6-22-76  pa 

Apilichlcoli,  Fli. 

0.00894 

0.161 

1.480 

0.160 

0.078 

0.114 

8-23-76  ia 

Apilichlcolc,  Fli. 

0.00894 

0.168 

1.627 

0.214 

0.042 

0.166 

9-12-76  mi 

Albuquerque,  N.M. 

0.00739 

0.024 

1.780 

0.111 

0.008 

0.100 

9-13-76  ia 

Albuquerque,  N.M. 

0.00740 

0.014 

0.334 

0.078 

0.021 

0.091 

9-12-76  la 

Albuquerque,  N.M. 

0.00744 

0.047 

-0.226 

0.107 

0.008 

0.076 

9-19-76  mi 

Albuquerque,  N.M. 

'0.00737 

0.027 

1.732 

0.107 

0.011 

0.127 

9-28-76  ai 

Odilon,  Utah 

0.00768 

0.033 

1.026 

0.061 

0.026 

0.077 

9-26-76  p> 

Ogden,  Utah 

0.00768 

0.033 

1.064 

0.047 

0.024 

0.080 

9-29-76  la 

Ogdai,  Utah 

0.00769 

0.032 

0.904 

0.041 

0.026 

0.069 

10-26-76  la 

Dayton,  Ohio 

0,00872 

0.07 6 

0.828 

0.036 

0.057 

0.112 

10-26-76  pa 

Dayton,  Ohio 

0.00872 

0.061 

1.074 

0.044 

0.053 

0.118 
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Equation  2-6 

Where:  K is  124.87  + 6.754  x 10-6  P 
s s 

P is  the  station  atmospheric  pressure  in  millibars 
(not  corrected  to  sea  level) 

b is  7.315  + . 000908P  + .0236T 

T is  the  station  air  temperature  in  degrees 
Celsius,  and 

Xj,  x„  are  the  elevations  above  sea  level  of  the 
* station  and  the  observer,  in  kilometers. 

Figure  2-7  shows  a typical  RAOB  density  profile  (see  x's)  and  the  Weibull 
function  fit.  A complete  table  of  station  pressure  observations  for  all  day- 
sites  of  the  project  can  be  found  in  Appendix  B of  Volume  2 of  this  report. 

The  total  vertical-beam-path  transmittance  due  to  pure  Rayleigh  scatter  is 
then  computed  as  follows: 

Equation  2-7 

2.4.2  Conversion  to  Sea-Level  Equivalent  Pressure 

As  stated  above,  the  SCAT3  model's  algorithms  were  derived  using  station 
pressure.  This  pressure  as  reported  by  weather  stations  is  usually  the 
equivalent  sea- level  pressure  obtained  by  reduction  of  the  station  pres- 
sure to  sea  level,  using  an  inverted  version  of  the  hypsometric  equation 
discussed  by  Gordon.9 

9 A.  H.  Gordon,  Elements  of  Dynamic  Meteorology,  D.  VanNostrand  Company, 

Inc,  New  York,  N.  Y.,  1962. 
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Figure  2-7.  Typical  Atmospheric  Density  Profile 
(23  Feb  76.  12002) 
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The  sea  level  pressure  is  the  input  parameter  for  SCAT3,  and  the  model 
converts  the  sea  level  pressure  to  station  pressure  by  means  of  the 
following  equation: 

14.  8251: 

Po  ■ P„  x 10  273.16  ♦ T Equation  2-8 

s 0 mv  ^ 


Where: 


Pq  is  the  sea  level  pressure  in  millibars 
Ps  is  the  station  pressure  in  millibars 
E is  the  station  elevation  in  kilometers,  and 
T is  the  mean  virtual  temperature  in  degrees  Celsius. 


In  a system  of  moist  air,  the  virtual  temperature,  Tmv>  is  the  temperature 
of  dry  air  having  the  same  density  and  pressure  as  the  moist  air  and  is  a 
function  of  station  pressure,  vapor  pressure,  and  temperature.  Because  in 
most  practical  cases,  the  virtual  temperature  is  not  significantly  dif- 
ferent from  the  actual  temperature,  they  were  considered  the  same  for  this 
study.  Another  approximation  that  was  made  was  that  the  station  elevation 
in  geopotential  meters  is  equal  to  the  elevation  in  meters. 


2.4.3  Volumetric  Scaling  for  Aerosol  Scattering 

Separation  of  the  measured  atmospheric  spectral  transmittance  into  its 
components  permitted  the  aerosol  extinction  coefficients  to  be  estimated 
for  each  day-site.  These  were  then  analyzed  as  functions  of  other  measurable 
parameters,  primarily  moisture  scale  height.  In  addition,  a joint  experiment 
was  conducted  in  Albuquerque,  N.M.,with  the  Atmospheric  Sciences  Laboratory  of 
the  U.S.  Army  Electronics  Command,  White  Sands,  N.M.,  in  which  measurements 
were  made  of  surface  dust  particle  concentrations.  The  experiment  was  con- 
ducted using  two  trailer  laboratories  located  together  in  an  area  of  sparse 
vegetation.  Aerosol  measurements  were  made  over  eight  consecutive  days 
(12-19  September  1976),  and  radiometric  readings  were  made  on  four  of  these 
days.  A description  of  the  experiment  and  a summary  of  the  data  is  contained 
in  para.  5.2,  Vol  2 of  this  report. 
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The  Albuquerque  experiment  was  used  to  refine  this  relationship  between 
the  aerosol  extinction  coefficients  and  moisture  scale  height  by  intro- 
ducing an  independent  parameter  for  dust  concentration  into  the  coefficient 
equations.  A complete  list  of  aerosol  extinction  coefficients  is  contained 
in  Table  2-1  above. 

The  equations,  as  modified  by  the  Albuquerque  experiment,  were: 

Wavelength  exponent,  C * 1.492-.000424Dq  ♦ .092MSH;  C > .5  Eq  2-9 

Extinction  Coefficient,  Kfl  - -.0705  ♦ .000015D0  ♦ .0831MSH;  Kfl  > .0481  Eq  2-10 

Where:  Dn  is  the  surface  dust  particle  concentration  in 

particles/cm3  equivalent  to  the  McClatchey  mode 
at  sea  level. 

0 < D0  < 13780,  if  1 < MSH  < 2 
Dq  « 1712,  if  MSH  < 1 or  MSH  > 2 

MSH  is  the  total  moisture  scale  height,  in  centimeters 
CO  < MSH  < 10) 

Figure  2-8  shows  the  results  of  these  equations  when  they  are  evaluated 
for  several  levels  of  moisture  scale  height  and  dust  concentration.  These 
parameters  are  used  to  estimate  the  aerosol  extinction  coefficient,  Ka,  and 
the  vertical  beam-path  transmittance  according  to  the  following  equation: 

Equation  2-11 

2.4.4  Altitude  Scaling  for  Aerosol  Scattering 

Altitude  scaling  of  the  aerosol  extinction  coefficients  was  possible  by 
using  the  assumption  that  the  amount  of  scattering  due  to  aerosols  would  be 
distributed  vertically  in  the  same  manner  as  the  amount  of  water  content. 

Thus,  the  extinction  coefficient  would  increase  approximately  exponentially 
toward  the  earth's  surface  and  approach  zero  in  space.  In  other  models, 
such  as  the  one  developed  by  McClatchey,  similar  distributional  shape  as- 
sumptions are  made, but  the  extinction  coefficient  distribution  is  not  asso- 
ciated with  that  of  moisture. 
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AEROSOL  EX 


l.S 


Note:  Numbers  given  at  data  points  are  precipitable  moisture  in  cm. 


Figure  2-8.  Aerosol  Coefficients  as  a Function  of  Dust  and 
Precipitable  Moisture 
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The  profiles  of  moisture  scale  height  versus  altitude  from  the  RAOB 
launches  were  used  to  model  the  vertical  distribution.  By  area-normalizing 
the  moisture-scale-height  curves  and  referencing  them  to  ground  level,  it 
was  possible  to  compute  a single-altitude  vector  that  would  reasonably  ap- 
proximate any  of  the  curves.  Figure  2-9  shows  three  typical  area-normal- 
ized distributions  that  have  a general  decreasing  exponential  shape  with  a 
pronounced  peak  at  about  9000  feet  above  ground  level  and  a lesser  peak 
at  about  22,000  feet.  Note  that  the  offset  in  the  peaks  in  Figure  2-9  is 
correlated  to  the  difference  in  altitude  above  sea  level  for  the  three  sites. 

By  scaling  the  moisture  scale  height  according  to  the  ratio  of  the  areas 
under  the  moisture  distribution  function,  the  scaled  values  of  the  aerosol 
extinction  coefficients  can  be  computed  using  the  same  equations  as  in 
para.  2.4.3.  Equation  2-11  could  then  be  used  to  compute  the  total  ver- 
tical beam-path  spectral  transmittance  due  to  aerosols. 


The  altitude  vector  for  the  area-normalized  moisture  scale  height  is  given 
in  Table  2-2.  Note  that  the  altitude  bands  are  not  in  even  increments. 

This  same  unequal  altitude  stratification  is  used  throughout  SCAT3  whenever 
altitude  scaling  is  estimated  from  discrete  data  rather  than  closed-form 
equations. 

2.4.5  Estimating  Moisture  Scale  Height  from  Surface  Meteorology 

The  on-board  weather  data  collection  system  and  the  NWS  observations  per- 
mitted an  examination  of  the  relationship  between  surface  synoptic  data  and 
detailed  atmpsoheric  profiles  from  RAOB's.  Because  it  was  anticipated  that 
the  total-column  moisture  scale  height  would  be  related  to  surface  absolute 
humidity,  a simple  mathematical  approximation  to  a psychometric  chart  was 
found  useful  in  converting  from  relative  humidity  and  temperature  to  absolute 
water  vapor  pressure. 


A nonlinear  equation  was  found  that  could  be  used  in  estimating  vapor  pres- 
sure from  temperature  (°F)  along  a fixed  relative  humidity  curve.  Relative 
humidity  therefore  enters  the  general  equation  for  vapor  pressure  as  a 
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TABLE  2-2 


PREC1PITABLE  MOISTURE  AS  A FUNCTION  OF  ALTITUDE 


ALT  (in  Kilometers) 


Area-Normalized  Mean  Row  Vector 


0 

- 

.3048 

.1588552 

.3048 

- 

.6096 

.1344122 

. 6096 

- 

.9144 

.1133070 

.9144 

- 

1.2192 

.9514909 

X 

1.2192 

- 

1.5240 

.8185085 

X 

1.5240 

- 

1.8288 

.6916665 

X 

1.8288 

- 

2. 1336 

.5781139 

X 

2.1336 

- 

2.4384 

.4907499 

X 

2.4384 

- 

3.0480 

.7593858 

X 

3.0480 

- 

3 . 6576 

. 5431665 

X 

3.6576 

- 

4.2672 

.3642236 

X 

4.2672 

- 

4.876S 

.2548728 

X 

4.8768 

- 

5.4864 

.1797324 

X 

5.48b4 

- 

6.0960 

.1170285 

X 

6.0960 

- 

7.6200 

. 1432864 

X 

7.6200 

- 

9.1440 

.4041271 

X 

9.1440 

- 

10.6680 

.1005052 

X 

10.6680 

- 

12.1920 

.2710833 

X 

12.1920 

- 

13.7160 

.1031096 

X 

13.7160 

- 

15.2400 

.6268420 

X 

15.2400 

- 

18.2880 

.8532016 

X 

18.2880 

- 

21.3360 

.5161841 

X 

21.3360 

- 

24.3840 

.3155043 

X 

24.3840 

- 

27.4320 

.1927280 

X 

27.4320 

- 

30.4800 

.1193245 

X 

multiplying  factor,  resulting  in: 

VP  = . 805RI!  x exp  f .1 


.183T-69  J 


Equation  2-12 


where  RU  is  the  relative  humidity,  expressed  as  a decimal,  and  T is  the 
temperature  in  degrees  F. 

In  the  above  equation,  the  vapor  pressure  is  given  in  inches  of  mercury, 
but  for  use  in  the  model,  the  equation  was  converted  to  vapor  pressure 
in  millibars  by  multiplying  by  33.9mb/inch  Hg.  In  addition,  the  input 
parameters  of  the  model  are  in  metric  units,  and  a temporary  conversion 
of  the  input  temperature  from  °C  to  °F  is  therefore  executed  internally 
in  the  SCAT3  program  to  make  it  compatible  with  Equation  2-12. 
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The  vapor  pressure  is  then  converted  to  precipitable  moisture  using  a linear 
relationship  that  exists  between  the  square  root  of  vapor  pressure  and  the 
natural  logarithm  of  the  precipitable  moisture.  This  relationship  had  been 
suggested  by  other  experimenters10’11,  and  also  appeared  to  fit  the  observa- 
tions made  on  this  project  although  the  exact  coefficients  deviate  from  the 
relationships  in  the  literature.  The  linear  fit  found  on  this  project  is 
shown  in  Figure  2-10  as  well  as  the  functions  suggested  in  the  two  references 
above. 

The  equation  used  to  compute  moisture  scale  height  (MSH)  is: 

MSH  = 0.134  exp  ^.659  • VP  5 J Equation  2-13 

where  VP  is  the  vapor  pressure  in  millibars. 

2.4.6  Scaling  of  Absorption  Transmittance 

Significant  absorption  bands  are  the  result  of  water,  ozone,  and  the  oxygen 
molecular  content  of  the  atmosphere.  The  resulting  absorption  log  trans- 
mittance can  be  expressed  in  the  following  vector  form,  as  was  discussed 
above  in  para.  2.3. 

logeTa(\)  = Kw  • VW(A)  ♦ Kq  • Vo  (A)  ♦ Ko  • VQ  (A)  Equation  2-14 

2 2 3 3 

Vectors  established  for  earlier  atmospheric  models  and  discussed  by  McClatchey12 
were  used  in  this  model  with  the  corresponding  coefficients,  Kw,  ko  , and  Kq 
determined  from  mobile  laboratory  transmittance  data  for  each  day-sfte. 

The  method  of  predicting  these  coefficients  is  discussed  below. 


10  J.L.  Montieth,  "An  Empirical  Method  of  Estimating  Long-Wave  Radiation 
Exchanges  in  the  British  Isles,"  Quarterly  Journal  of  the  Royal  Meteo- 
rological Society,  551.521.32. 

11  S.B.  Idso,  "Atmospheric  Attenuation  of  Solar  Radiation,"  Journal  of 
Atmospheric  Sciences,  Vol.  26,  pp  1088-1095. 

12  R.A.  McClatchey  et  al.  "Optical  Properties  of  the  Atmosphere,"  Air  Force 
Cambridge  Research  Laboratories  -72-0497,  Aug  72,  Environmental  Research 
Papers  »411. 
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2.4.6. 1 Water  Vapor  Absorption 


Water  absorption  occurs  in  four  distinct  wavelength  bands:  720-735,  810-840, 
895-990,  and  1105-1150  nanometers.  The  depth  of  each  band  is  given  by  the 
product  of  the  water  absorption  vector,  Vw(A),  and  a scaling  coefficient  »kw, 
found  to  be  related  to  moisture  scale  height.  The  water  absorption  vector  is 
shown  in  Figure  2-11.  This  vector  was  fitted  to  the  measured  absorption  log 
transmittance  data,  using  the  method  of  least  squares,  to  determine  the  value 
of  Kw  for  each  day-site.  Because  instrument  non-linearities  were  suspected 
beyond  1 pm,  only  the  first  three  bands  were  used  in  those  regressions. 

As  was  anticipated,  the  relationship  between  the  vector  coefficient,  K^,  and 
the  moisture  scale  height  was  not  linear.  Figure  2-12  shows  the  regression 
equation  fitted  to  the  estimated  coefficients  as  a function  of  moisture  scale 
height.  This  regression  equation  is: 


1^  * -0.0645  MSB  0,713  Equation  2-15 

The  model  also  scales  the  water  bands  as  a function  of  the  observer's  altitude 
and  the  terrain  altitude.  This  is  done  by  ratioing  the  appropriate  partial 
area,  from  the  ground  elevation  to  the  observer's  altitude,  to  the  total  area 
under  the  moisture  scale  height  as  a function  of  altitude  (see  Figure  2-9  for 
examples) . 


2.4.6. 2 Oxygen  Absorption  Line 


A molecular  oxygen  absorption  line  has  been  modeled  in  the  760-765  nm  band. 

As  was  the  case  with  other  absorption  transmittance  data,  this  narrow  absorp- 
tion band  was  fitted  using  an  arbitrary  vector,  V (A)  and  coefficient,  V . 

°2  2 

Analysis  of  all  day-site  data  resulted  in  an  average  value  of  0.0905  for 

K and  a standard  deviation  of  0.042.  Since  the  absolute  amount  of  oxygen 
©2 

is^relatively  constant,  it  was  presumed  that  the  variance  about  the  mean  Kq 
was  the  result  of  pure  measurement  error;  it  was  therefore  not  modeled  as 


RELATIVE  ABSORPTANCE 


X 


0 .50  1.00  1.50  2.00  2.50  3.00  3.50  4.00  4.50  5.00  5.50 

SCALE  HEIGHT  (in  cm) 


.Figure  2-12.  Model  Equation  Fitted  to  Water  Absorption 
Coefficient  Versus  Moisture  Scale  Height 
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I 


a function  of  any  other  measurable  parameter.  The  vector  has  the  following 
values  at  the  wavelengths  listed  below. 


350-755 

nm 

.00 

760 

nm 

1.00 

765 

nm 

.24 

770-1150 

nm 

.00 

Altitude  scaling  can  be  accomplished  for  the  oxygen  line  if  one  assumes  that 
the  absorption  is  linear  with  atmospheric  density.  Thus,  the  absorption  for 
a partial  atmosphere  can  be  computed  by  integrating  the  atmospheric  density 
profiles  in  the  same  manner  as  for  molecular  scattering  discussed  above  in 
para.  2.4.1. 


2. 4. 6. 3 Ozone  Absorption 


Molecular  ozone  absorption  occurs  in  a broad  band  in  the  visible  portion  of 
the  spectrum  between  470  and  710  nm.  The  amount  of  this  absorption  could  be 
estimated  from  analysis  of  actual  path  transmittance  measurements  (see  para. 

2.3).  These  measured  absorption  spectra  were  fitted  with  the  vector  Vo  (X) 
and  the  coefficient  K by  regression  techniques.  The  ozon6  vector  is  ^ 
shown  in  Figure  2-13.  The  average  of  the  ozone  coefficients  was  -0.049 
and  the  standard  deviation  was  0.044. 


Unlike  oxygen  absorption,  which  would  be  expected  to  vary  only  slightly  with 
physical  changes  in  the  atmosphere,  large  excursions  occur  seasonally  and  by 
latitude  in  ozone  concentration.  In  a previous  study  by  Craig13,  the  ozone 
scale  height  (in  mm)  as  distributed  seasonally  in  the  northern  hemisphere  was 
used  to  scale  the  absorption  transmittance.  The  absorption  was  assumed  to  be 
a gaussian  function  having  a peak  value  of  0.95  at  590  nm  for  a scale  height  of 
of  3.8  nm.  When  logged  and  peak-normalized  to  1.0,  this  function  became  the 
ozone  absorption  vector  used  in  this  study  as  can  be  seen  in  Figure  2-13. 

The  scaling  of  ozone  absorption  can  also  be  accomplished  as  a function  of 
altitude  according  to  the  ozone  concentration  distribution  from  the  Handbook 


13  R.A.  Craig,  "The  Upper  Atmosphere,  Meteorology,  and  Physics,  " Academic 
Press,  1965,  pp  179-181. 
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of  Geophysics  and  Space  Environments.  This  distribution  is  shown  in 
Figure  2-14.  The  integral  of  the  distribution  function  between  the  obser- 
ver's altitude  and  the  terrain  altitude  ratioed  to  the  total  earth- to-space 
scale  height  (3.2mm)  results  in  the  scaled  concentration  from  which  the 
vector  coefficient  can  be  computed. 

A proportionality  constant  CQ  existed  between  the  vector  coefficient  VQ  and 

the  ozone  scale  height,  as  reported  by  Craig.13  This  proportionality  constant 
can  be  expressed  as  follows: 


l°g« 


1-.05  Exp 


(X- .59) 2 
.00269 


(X) 

T59) 


Equation  2-16 


Evaluating  equation  2-16  at  the  normalizing  wavelength  .59  nm  results  in  the 

proportionality  constant  CQ  , which  equals  -.0513. 

3 


If  the  ozone  coefficient  is  scaled  linearly  with  the  ozone  scale  height 
relative  to  the  base  level  of  3.8mm  and  the  constant  of  proportionality  is 
introduced,  the  following  equation  results. 


K *=  C • * - -0135  OZ  Equation  2-17 

°3  °3  3‘8 

where  OZ  is  the  ozone  scale  height  in  mm,  which  is  consistent  with  the  data 
and  definition  discussed  by  Craig.13  Scaling  in  the  SCAT3  model  is  accom- 
plished using  the  reference  distribution  of  ozone  concentration  as  a function 
of  latitude  and  date.  The  actual  distribution  of  the  coefficient  K as  ob- 

°3 

served  during  this  project  is  shown  in  Figure  2-15. 
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OZONE  COEFFICIENT  (-K 
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Figure  2-15.  Distribution  of  Ozone  by  Date 
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2.4.7  Scaling  for  Observer's  Zenith  Angle 


Atmospheric  transmittance,  of  both  scattering  and  absorption  types,  is 
assumed  to  vary  with  the  observer's  zenith  angle  according  to  the  secant  law 
as  discussed  by  Robinson14  except  at  low  solar  altitudes  (below  30°)  where 
the  air  mass  numbers  of  Bemporad15  are  used.  Analysis  of  the  spectral  trans- 
mittance data  indicated  that  the  total  extinction  coefficient  was  linear  with 
air  mass  when  computed  in  this  manner.  Water  absorption  level  data  tended  to 
be  more  erratic  in  this  relationship,  probably  due  to  diurnal  variations  in 
atmospheric  moisture  profiles;  however,  for  the  most  part  even  these  bands 
were  well  approximated  by  linear  air  mass  scaling.  SCAT3  uses  the  following 
equation  for  this  purpose: 


te(x)  = tvcx)m 


Where:  t (X)  is  the  total  vertical  spectral  transmittance  of  the 

atmosphere  (after  scaling  for  altitude,  moisture  scale 
height,  and  dust). 

tg (X)  is  the  total  spectral  transmittance  of  the  atmosphere 
for  the  angle  0. 


6 

m 


is  the  observer's  zenith  angle,  and 


is  the  air  mass 


(sec  0 for  the  range  0°  < 0 < 60°, 

Bemporad  numbers  for  the  range  0 > 60°. 


Table  2-3  lists  the  Bemporad  air-mass  numbers15  as  a function  of  the  observer's 
zenith  angle  in  1-degree  increments. 


14  N.  Robinson,  Solar  Radiation,  Solar  Physics  Laboratory,  Israel  Institute 
of  Technology,  Elsevier  Publishing  Co.,  1966,  p 52. 

15  A. Bemporad,  "Search  for  a New  Empirical  Formula  for  the  Representation 

of  the  Intensity  of  Solar  Radiation  with  Zenith  Angle,"  Meteorologische 
Zeitschrift,  Vol  24,  July, 1907.  ~~ 
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TABLE  2-3 

BEMPORAD  AIR-MASS  NUMBERS 
(Observer's  Zenith  Angle  > 60®) 


Zenith  Angle 
(in  degrees) 

Air  Mass 

Zenith  Angle 
(in  degrees) 

Air  Mass 

61 

2.06 

76° 

4.07 

62 

2.12 

77 

4.37 

63 

2.19 

78 

4.72 

64 

2.27 

79 

5.12 

65 

2.36 

80 

5.60 

66 

2.45 

81 

6.18 

67 

2.55 

82 

6.88 

68 

2.65 

83 

7.77 

69 

2.77 

84 

8.90 

70 

2.90 

85 

10.39 

71 

3.05 

86 

12.44 

72 

3.21 

87 

15.36 

73 

3.39 

88 

19.79 

74 

3.59 

89 

26.96 

75 

3.82 

90 

39.65 

2.5  THE  RELATIONSHIP  BETWEEN  METEOROLOGICAL/CLIMATOLOGICAL  DATA  AND  T550 

As  stated  above,  the  atmospheric  model  was  designed  to  use  temperature,  pres- 
sure, and  relative  humidity  as  basic  input  parameters.  Situations  may  occur, 
however,  where  it  would  be  desirable  not  to  use  the  surface  meteorology  data 
and  key  the  model  to  a given  transmittance  value.  For  this  purpose,  the  tran- 
mittance  at  550  nm  was  arbitrarily  selected  as  the  alternative  input  parameter. 
Provision  was  therefore  made  in  the  SCAT3  model  to  input  a T550  value  instead 
of  the  meteorological  parameters.  When  the  T550  input  option  is  used,  the 
model  creates  a set  of  meteorological  parameters  that  will  give  the  specified 
transmittance  value  at  550  nanometers.  When  the  created  values  of  temperature, 
pressure,  and  relative  humidity,  are  used,  the  model  continues  its  calculation 
of  the  spectral  transmittance  in  the  usual  manner. 


To  create  the  meteorological  parameters,  the  model  assumes  a standard 
pressure  of  1013.25  millibars  and  a standard  temperature  of  20°C  unless 
they  are  specified  otherwise  by  the  program  user.  SCAT3  then  solves  the 
equations  in  reverse  order  to  calculate  the  corresponding  relative  humidity. 
This  is  done  by  first  calculating  the  Rayleigh  coefficient,  Kr,  using  the 
following  equation: 


.008925 

1013.25 


x AP 


Equation  2-18 


In  this  calculation,  the  aerosol  coefficient  C is  taken  to  be  the  overall 
project  average  value  of  1.1014  and  if  the  values  of  Kj.,  C,  t are  combined, 
Equation  2-3  can  be  solved  for  the  aerosol  coefficient  K.a  after  dividing  out 
an  assumed  amount  of  ozone  absorption  transmittance  that  is  derived  from  an 
input  or  default  value  for  ozone  scale  height. 

To  compute  the  equivalent  aerosol  scattering  coefficient,  the  following  equa- 
tion is  used: 


K * 0.545  • lnT550  ♦ 5.960  K_  ♦ .296  K 

8 * C 


Equation  2-19 


Given  K and  C,  the  moisture  scale  height  can  be  computed  using  the  inversion 
8 

of  equation  2-10,  as  follows: 


MSH  ■ 11.57  K ♦ . 420C  ♦ .189 

a 


Equation  2-20 


The  vapor  pressure  in  millibars  can  then  be  computed  using  the  inversion  of 
equation  2-13,  as  shown  below. 


VP  , ( In (MSH) ♦ 2.011  V 
\ .659  / 


Equation  2-21 


As  a final  step,  the  relative  humidity  is  computed  using  the  assumed  temper- 
ature of  20°C  and  VP  in  inverted  equation  2-12. 


RH  - 1.242  VP  exp  -.183T’ 


Equation  2-22 


I 


If  the  relative  humidity  when  computed  is  greater  than  1.0,  then  the 
program  will  increase  T in  one-degree  increments  until  a relative  humidity 
of  less  than  1.0  is  achieved.  The  algorithm  then  satisfies  the  require- 
ment of  T550  while  physically  possible  surface  meteorological  conditions 
are  selected. 

As  can  be  seen  in  the  description  above,  when  T550  is  specified,  average 
values  of  the  two  meteorological  parameters  temperature  and  pressure  are 
assumed  when  the  third,  relative  humidity,  is  calculated.  This  approach  is 
based  on  the  assumption  that  relative  humidity  changes  through  a wide 
range  of  values  in  a way  that  is  harder  to  approximate  accurately  than 
by  using  values  of  temperature  or  pressure.  In  addition,  pressure  variations 
affect  the  transmittance  less  than  changes  in  temperature  and  relative  humid- 
ity do  (see  Figures  D-l  to  D-3  in  Appendix  D) . If  it  is  known  that  for  a 
particular  case  the  temperature  will  be  significantly  different  from  20°C, 
the  user  may  choose  to  input  T. 


For  comparison  purposes,  measured  temperature,  pressure,  and  relative  humidity 
were  used  along  with  transmittance  values  at  550  nm  to  reconstruct  the  spec- 
tral transmittances  in  Figures  2-i6  and  2-17.  Figure  2-16  shows  the  compar- 
ison for  a clear  day  having  a low  moisture  scale  height;  Figure  2-17,  on  the 
other  hand,  shows  a similar  comparison  for  a high-moisture,  hazy  day.  Figure 
2-16  shows  practically  no  difference  in  the  transmittance  generated  by  either 
method  for  low-moisture  days.  A slightly  larger  but  still  small  difference 
in  transmittance  across  the  spectrum  is  evident  in  Figure  2-17  for  high- 
moisture  days.  Because  these  differences  are  small,  it  was  concluded  that 
the  model  can  be  used  successfully  in  either  the  T550  or  the  meteorological 
mode. 


To  make  the  relationship  between  relative  humidity,  temperature,  and  T550 
easier  to  understand.  Figure  2-18  was  prepared  by  operating  the  program  over 
a large  range  of  parameter  values.  In  the  figure,  lines  of  constant  T550  were 
interpolated  to  show  the  trade-offs  in  selection  of  temperature  and  relative 
humidity. 
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Figure  2-16.  Comparison  of  T550  and  Meterologically  Generated  Spectral 
Log  Transmittance  for  a Low-Moisture  Day 
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Figure  2-17.  Comparison  of  T550  and  Meterologically  Generated  Spectral 
Log  Transmittance  for  a High-Moisture  Day 
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SECTION  3 


ANALYSIS  AND  MODELING  OF  THE  ATMOSPHERIC  PATH  RADIANCE 


I 

[I 

r 


3.1  NORMALIZATION  OF  MEASURED  DATA 


The  atmospheric  path  spectral  radiance  measured  by  the  sky  spectroradiometer 
was  analyzed  in  the  sequence  of  manual  and  computer  operations  referred  to  in 
Section  1.  Each  individual  spectral  measurement  was  normalized  to  a common 
path- length  and  the  geometry  that  corresponds  to  an  observer  looking  vertically 
downward  from  outside  the  atmosphere.*  This  was  accomplished  as  a two-stage 
correction,  the  first  compensating  for  total  optical  path  length,  and  the 
second  correcting  for  the  relative  fraction  scattered  as  a function  of  obser- 
ver/ sun  geometry.  The  measured  spectral  path  radiance  was: 

N ( A , , 4>^ » 0)  = Measured  sky  radiance  in  watts/m2 
ster/5  nm 


Where:  X 
* 

*b 

e 


is  the  wavelength  in  nanometers 
is  the  observer's  zenith  angle 

is  the  angle  between  the  observer's  ray  and  the  sun's 
ray  measured  relative  to  pure  backscatter,  and 

is  the  solar  zenith  angle  (90°  minus  solar  altitude) . 


The  sky  spectroradiometer  telescope  was  always  pointed  at  a specific  point 
in  the  sky,  the  azimuth  (AZ)  and  elevation  (EL)  of  which  maintained  a pre- 
cise relationship  to  the  sun's  position  as  discussed  by  Edgerton.16  For 
positions  above  and  below  35°  solar  altitude,  the  relationship  was  different 
and  resulted  in  the  following  equations: 


Measurement  angle  $ 


18O°-20;  55°  < 0 < 90° 
9O°-0;  0 < 55° 


Equation  3-1 


* Another  correction  was  required  at  extreme  low  solar  altitudes  to  com- 
pensate for  the  non- instantaneous  nature  of  the  spectral  radiometric 
measurements  referred  to  in  Appendix  E. 

16  C.F.  Edgerton,  Relationship  Between  Meteorological  Conditions  and  Optical 
Properties  of  the  Atmosphere,  U.  of  Cal.  Scripps  Inst,  of  Oceonography 
and  Visibility  Lab,  May  1967. 


a 
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Measurement  angle 


6 

90* 


55°  < 6 < 90° i 
e < 55° 


Equation  3-2 


Measurement  azimuth  A2  = 180°  - Solar  Azimuth 


The  above  restrictions  in  collection  geometry  significantly  simplified 
normalization  of  the  spectral  sky  radiance  data  by  eliminating  the  re- 
quirement for  scatter-angle  correction  for  most  of  the  measurements. 
Figure  3-1  shows  the  geometry  used  to  normalize  measurements  to  the  verti- 
cal look. 


Figure  3-1.  Measurement  of  Sky  Radiance 
(Sun's  and  Observer's  Rays  are  Coplanar) 


f 


The  two  equations  of  normalization  are  shown  below.  For  solar  zenith  angles 
between  55°  and  90°: 


N(A,O°,0,0) 


N(A,  18O°-20,0,0) 


/ i-yx) 

\ l-t"(A) 


) 


Equation  3-3 


Where:  ts(A)  is  the  vertical  atmospheric  path  transmittance 
due  to  scattering,  and 

m is  the  air  mass  corresponding  to  a measurement 
zenith  angle  of  $ = 180-20 


For  solar  zenith  angles  less  than  55°: 

1-t  fA)  fm 

N(A,O°,0,0)  = N(A,9O°-0,9O°,0)  • - . Equation  3-4 

l-tm(A)  f(90  ' 

s 

Where:  m is  the  air  mass  corresponding  to  a measurement 
zenith  angle  of  $ = 90-8,  and 

f (0)  ,f  (90°)  are  the  polar  scatter  function  evaluated  at  4>,  = 8 
and  4>  = 90°,  respectively.* 

As  a result  of  the  uniformity  restrictions  imposed  by  Equation  3-4  assump- 
tions, not  all  the  spectral  radiance  data  could  be  used  in  this  study. 

Although  data  from  the  days  having  obvious  cloud  structure  were  avoided,  a 
satisfactory  range  of  meteorological  conditions  were  included.  In  all,  20 
of  the  day-site  measurement  series  over  a range  of  solar  altitudes  were 
examined  that  had  moisture  scale  heights  from  0.23  to  4.2  cm  and  terrain  re- 
flectances from  forest  to  snow  cover. 

3.2  Vertical  Path  Radiance  Analysis 

As  is  shown  in  Equations  3-3  and  3-4,  once  the  path  radiance  was  normalized, 
it  was  strictly  a function  of  wavelength  and  solar  zenith  angle,  0.  Normal- 
ization of  these  data  in  turn  by  the  extraterrestrial  irradiance  and  vertical 


* For  the  derivation  of  the  polar  scatter  function,  see  Section  5 of 
this  report. 


3-3 


path  absorption  transmittance  permitted  the  examination  of  vertical  path 
radiance  as  a function  of  the  total  fraction  scattered,  as  estimated  by 


1 - ts(X).  The  mathematical  model  below  for  path  radiance  was  hypothesized 


and  then  regressed  against  all  available  spectral  vertical  path  radiance 
data. 

b2 

N(X  ,6)  ■ b,Ho/%(X)  | l-te (X)  | t (X)  Equation  3-5 


>lHsc<X>  1-ts(X) 


Where:  H (X)  is  the  spectral  extraterrestrial  irradiance 
in  watts/m2. 


t (X) • t (X)  is  the  vertical  path  transmittance  due  to 
scattering  and  absorption,  and 


bj.b^  are  regression-determined  parameters. 


Dividing  Equations  3-5  through  by  the  solar  constant,  Hs(.(A),  and  the 


absorption  transmittance,  t,(A),  creates  a normalized  radiance  variable 

Cl 


as  a function  of  soatter  fraction. 


N(x,e) 

"sc'1) 


b,  [1  - t5w] 


Equation  3-6 


3.2.1  The  Shape  Parameter,  b2 


Through  study  of  the  log  of  Equation  3-6,  the  relationship  of  the  shape 
parameter  b2  to  solar  altitude  and/or  moisture  scale  height  was  determined. 

A typical  example  of  this  relationship  is  solar  altitude  shown  in  Figures 
3-2  and  3-3.  On  the  basis  of  examination  of  all  data,  a two-piece  linear 
model  was  fitted  to  these  data  as  shown  in  Figures  3-2  and  3-3.  In  addition, 
a slight  correlation  was  apparent  with  moisture  scale  height  as  can  be 
seen  in  Figures  3-4  and  3-5,  and  the  equations  were  modified  to  incorporate 
this  minor  effect.  The  resulting  equation  used  to  compute  the  shape  param- 
eter was: 


b_  (shape  parameter)  * 1.037  ♦ .015  SA  ♦ .0030  SA 
2 (MSH  - 1.0) 


Equation  3-7 


Where:  SA  ■ Solar  Altitude  (#) - 


SA,  SA  < 40* 
40* , SA  > 40* 


and 


ii 
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SOLAR  ALTITUDE  (in  degrees) 


Figure  3-2.  B2  Coefficient  as  a Function  of  Solar  Altitude 
for  June  11,  1976  Data 


SOLAR  ALTITUDE  (in  degrees) 

Figure  3-3.  B2  Coefficient  as  a Function  of  Solar  Altitude 
for  August  li,  1976  Data 
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MOISTURE  SCALE  HEIGHT  IN  CBJTIHETERS 


Figure  3-4.  B2  Coefficient  As  a Function  of  Moisture  Scale  Height 
for  All  Day-Site  Data  (SA  * 20°) 


MOISTURE  SCALE  HE! (BIT  IN  CHITIMETERS 


Figure  3-5.  B2  Coefficient  As  a Function  of  Moisture  Scale  Height 
for  All  Day-Site  Data  (SA  - 40°) 

3-6 


1.0,  MSH  < l.o 


MSH  * Moisture  Scale  Height  (cm)  «<  MSH,  1.0  s MSH  i 4.0  ) 

I 4.0,  MSH  > 4.0  I 

As  Equation  3-7  indicates,  the  shape  coefficient  is  only  allowed  to 
vary  for  solar  altitudes  below  40°  and  for  moisture  scale  heights  between 
1 and  4 cm.  These  are  empirical  limits  imposed  from  observations  made  from 
plots  of  b2  versus  solar  altitude  and  moisture  scale  height. 

3.2.2  The  Magnitude  Parameter,  b^ 


Once  a suitable  algorithm  was  found  for  the  shape  coefficient,  b2,  it  was 
possible  to  investigate  the  magnitude  coefficient,  bj.  Linear  regressions 


were  run  against  the  normalized  measured  spectral  atmospheric  radiance  data 


using  Equation  3-6.  This  time,  however,  the  equation  was  not  logged  and 
the  values  of  b2  were  fixed  as  calculated  by  Equation  3-7.  The  series  of 
scalars,  bj,  developed  for  the  normalized  data  for  each  day  were  then 
studied  as  functions  of  solar  altitude,  moisture  scale  height,  and  terrain 
reflectance.  No  spectral  shaping  could  be  observed  in  the  normalized  path 
radiance  data  that  could  be  correlated  to  terrain  spectral  reflectance. 


The  terrain  reflectance  at  650  nm;  that  is,  R^q. 
single  descriptor  for  correlation  studies  with  bj 


was  therefore  used  as  a 


The  equation  given  below  was  fitted  to  plots  of  bj  versus  solar  altitude: 


bx  * Cj  sin  SA  3°  < SA  < 90°  Equation  3-8 

The  coefficient  Cj  was  found  to  correlate  with  moisture  scale  height,  MSH, 
and  the  terrain  reflectance  descriptor,  R650»  as  follows: 

Cj  « .2158  ♦ .0978R65q  - .0298MSH  Equation  3-9 

Figures  3-6  through  3-8  show  three  typical  examples  of  the  bj  coefficient 
plotted  versus  solar  altitude. 


■■■—  w «i  mnrnmmmmmmmmmtitmm 


(MSH*.47  CM,  Snow  Surround  R=.60) 


Equational  Fit 


SOLAR  ALTITUDE 

Figure  3-6.  Dependence  of  the  B.  Coefficient  for  Feb  24,  1976 


(MSH= . 76CM,  Snow  Surround  R«.20) 


Equational  Fit 


SOLAR  ALTITUDE 


Figure  3-7.  Dependence  of  the  B.  Coefficient  for  Jan  11,  1976 


(MSH=2. 91  CM,  Forest  Surround  R=.08) 


Equational  Fit 


SOLAR  ALTITUDE 


Figure  3-8.  Dependence  of  the  B,  Coefficient  for  Aug  11,  1976 


3.2.3  Very  Low  Solar  Altitude  Modeling 


Only  a few  sky  radiance  measurements  were  available  below  three  degrees 
solar  altitude.  To  continue  the  model's  applicability  down  to  minus  five 
degrees  solar  altitude,  it  was  necessary  to  adopt  a relative  scaling  poly- 
nomial for  sky  radiance  that  was  independent  of  moisture  scale  height  or 
surround  reflectance.  The  polynomial  itself  does  not  involve  spectral 
shaping;  in  addition,  it  contains  no  parameters  other  than  solar  altitude 
and  is  used  strictly  as  a relative  magnitude  adjustment  for  the  sky  radi- 
ance term.  The  adjustment  polynomial  P(SA)*  had  been  evaluated  at  three 
degrees  solar  altitude,  and  this  value  was  divided  into  the  coefficients 
to  normalize  the  polynomial  to  the  lowest  valid  calculated  spectral  sky 
radiance.  Equation  3-5  for  the  model  is  evaluated  for  the  appropriate 
condition  of  moisture  scale  height  and  terrain  reflectance  at  three  de- 
grees solar  altitude;  the  spectral  sky  radiance  is  then  multiplied  wave- 
length-by-wavelength by  the  normalized  adjustment  polynomial  evaluated  at 
the  solar  altitude  of  interest.  The  polynomial  and  scaling  approach  are 
explained  in  Appendix  F.  The  normalized  factor  used  for  spectral  sky 
radiance  is  0.01975. 

Because  refraction  of  the  sun's  rays  in  the  atmosphere  becomes  substantial 
below  five  degrees  solar  altitude,  a correction  must  be  applied  to  solar 
altitude  as  used  in  Equation  3-8.  The  corrected  solar  altitude  is  re- 
ferred to  as  the  "apparent"  solar  altitude  and  is  described  mathematical ly 
in  Appendix  G.** 


* See  para.  4.4  and  Appendix  F for  a description  of  the  polynomial 
P(SA) . 

*w  See  Appendix  G for  the  correction  applied  to  solar  altitude  below  5° 
to  compensate  for  refraction  of  the  sun's  rays  in  the  atmosphere. 
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3.3  ATMOSPHERIC  PATH  RADIANCE  MODEL  IN  SCAT3 


Starting  with  the  three  parameters  solar  altitude,  moisture  scale  height, 
and  the  terrain  reflectance  at  650  nm,  the  vertical  atmospheric  path  radiance 
can  be  computed  as  follows  from  Equation  3-5  if  the  scattering  and  absorp- 
tion transmittance  ts(A)  and  ta(A)  are  first  computed  using  the  equations  dis- 
cussed in  Section  2 of  this  report. 

b2 

N(A,0°,e,6)  = b.H  (A)  [1  - t (A)  ] t(A)  Equation  3-10 

1 SC  S 3 

Where  6 is  90°  minus  the  solar  altitude. 


Scaling  for  the  path  length  and  relative  angle  from  back  scatter  results 
in  the  following  general  expression  for  any  path  through  the  atmosphere: 


N(A,*,*b,6)  = Ks*  bx (e»MSH,R65Q)  . Hsc(A)  • [1  - ts(A)] 


b2(0,MSH) 


VA) 


1 - VX) 


Equation  3-11 


Where: 


A is  the  wavelength  in  nanometers 

6 is  the  solar  zenith  angle 

is  the  observer's  zenith  angle 

4>.  is  the  relative  angle  between  direct  backscatter 
and  the  observer's  line-of-sight 

MSH  is  the  moisture  scale  height  in  cm 


M($)  is  air  mass  corresponding  to  the  observer's 
zenith  angle 

f ($, A)  is  the  polar  scatter  function  evaluated  at  the 
observer's  zenith  angle 


^(♦u.A)  is  the  polar  scatter  function  evaluated  at  the  relative 
angle  from  backscatter,  and 

Ks  is  1.0  when  SA>3#,  or  is  p(SA)  when  SK<3*. 
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3.4  ANALYSIS  OF  THE  VOLUMETRIC  SCATTERING  FUNCTION,  f(*,X)* 


3.4.1  Data  Base 

An  estimate  of  the  volumetric  scattering  function  (the  relative  intensity 
of  scattered  light  as  a function  of  the  angle  around  a volume  of  atmo- 
sphere illuminated  by  a collimated  source)  was  sought  that  could  be  made 
to  vary  spectrally  and  according  to  atmospheric  optical  conditions.  For 
this  purpose,  sky  spectral  radiance  measurements  were  made  during  the  data 
collection  program  in  15-degree  increments  of  azimuth  between  0 and  360 
degrees,  and  at  an  elevation  of  7 degrees  above  the  horizon**.  Although  the 
complete  hemisphere  of  downward  radiance  was  actually  mapped  at  15-degree 
increments  of  elevation  from  7 degrees  to  83  degrees  above  the  horizon,  only 
the  7 degree  ring  was  used  in  this  analysis  because  this  elevation  ring 
measured  nearest  the  true  profile  of  the  scattering  function  when  the  sun  was 
near  the  horizon. The  data  base  was  also  culled  by  eliminating  all  measure- 
ments made  above  10  degrees  solar  altitude.  The  spectral  dependence  was 
evaluated  by  making  measurements  at  the  following  three  wavelengths:  450 
nanometers,  650  nanometers,  and  850  nanometers. 

Because  the  measurements  were  made  over  a finite  time  and  from  the  ground 
looking  out  toward  space,  corrections  had  to  be  applied  to  the  measurements 
to  convert  them  into  equivalent,  instantaneous,  gonioradiometric  measurements 
of  an  equivalent  volume  of  atmosphere.  The  method  of  correcting  the  data 
and  its  subsequent  reduction  to  equational  form  is  the  subject  of  this  section 
and  Appendix  H of  this  volume. 

3.4.2  Time  Correction 

At  the  same  time  data  was  being  collected  by  the  sky  radiometer  at  different 
instrument  azimuths  and  elevations,  data  was  collected  for  horizontal  sky 
irradiance  by  one  of  the  two  irradiometers.  Since  this  irradiometer  was 
oriented  in  the  same  position  for  all  readings  and  at  the  same  wavelength  as 


* See  para.  3.3  for  explanation  of  the  use  of  this  function. 

**  See  Section  3,  Vol  2 for  an  explanation  of  the  polar  scattering 


measurements . 


J 


the  sky  radiometer,  any  changes  in  this  measurement  were  changes  corres- 
ponding to  changes  in  sky  conditions.  The  simultaneous  sky  irradiance  data, 
taken  at  a specific  azimuth,  were  divided  into  the  sky  irradiance  data  cor- 
responding to  the  time  of  the  sky  radiometer  direct  backscatter  measurements 
made  at  zero  degrees  azimuth  from  the  sun.  This  fraction  was  then  multiplied 
by  the  sky  radiance  at  the  corresponding  azimuth  to  give  the  following  time- 
corrected  radiance. 


n' (AZ) 


N(AZ) 


H (0) 
H(AZ) 


Where:  N (AZ)  is  the  time-corrected  sky  radiance  at  azimuth 
angle  AZ 

N(AZ)  is  the  sky  radiance  at  angle  AZ  before  time  cor- 
rection. 


H(0)  is  the  sky  irradiance  data  recorded  at  the  same 
time  as  the  sky  radiance  for  pure  backscatter, 
made  at  zero  degrees  azimuth  from  the  sun. 

H(AZ)  is  the  sky  irradiance  data  recorded  at  the  same 
time  as  the  sky  radiance  at  angle  AZ  f and 

AZ  is  the  instrument  azimuth  angle  measured  relative 

to  the  solar  azimuth  (zero  degrees  is  directly 
away  from  sun) . 


3.4.3  Path  Length  Correction 


As  stated  above,  only  solar  altitude  data  less  than  10  degrees  and  at  a 7- 
degree  instrument  elevation  were  chosen  for  use.  In  the  15-degree  azimuth 
increments,  there  are  24  points  in  the  7 degree  ring,  each  of  which  is  illum- 
inated differently,  due  to  its  atmospheric  path  length  relative  to  the  sun. 
Since  the  atmosphere  absorbs  and  scatters  a fraction  of  the  incident  light, 
the  longer  the  path  the  light  must  travel,  the  smaller  the  amount  of  light 
that  reaches  a particular  point.  The  amount  of  light  reaching  its  destina- 
tion depends  on  the  length  of  the  path  that  the  light  must  travel  and  the 
extinction  coefficient  that  determines  the  amount  of  light  the  atmosphere  ab- 
sorbs or  scatters. 
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Similarly,  the  total  energy  reaching  the  sky  radiometer  is  the  integration 
of  scattering  gains  and  losses  along  the  path  toward  the  radiometer  and 
along  its  line  of  sight.  The  double  integration,  sun-to-volume  and  volume- 
to-radiometer,  was  digitally  approximated  using  a Hewlett-Packard  9830 
program,  the  arguments  for  which  were  the  measurement's  starting  solar  geom- 
etry and  the  vertical  air-path  transmittance  as  determined  from  the  day- 
site  measurements.*  In  the  program,  an  azimuth-angle-dependent  correction 
was  computed  that  was  applied  to  the  sky  radiometer  measurements.  Details 
about  the  deviation  of  this  correction  are  discussed  in  Appendix  H of  this 
volume. 

Correction  of  the  path  length  was  made  relative  to  the  solar  azimuth,  as 
was  the  time  correction.  Twenty- four  correction  values  were  computed  and 

divided  into  the  value  for  the  azimuth  that  corresponded  to  backscatter. 

These  in  turn  were  multiplied  by  the  data  already  corrected  for  time  as  indi- 
cated in  the  equation  above,  resulting  in  the  following: 

N"(AZ)  * N ' (AZ)  • Equation  3-12 


Where:  N"(AZ) 

is  the  path- length  and  time- corrected  sky  radiance 
at  azimuth  angle  AZ 

N' (AZ) 

is  the  sky  radiance  at  azimuth  angle  AZ  before 
path  length  correction 

PL(0) 

is  the  relative  amount  of  light  reaching  the 
observer  at  the  azimuth,  corresponding  to  back- 
scatter 

PL(AZ) 

is  the  relative  amount  of  light  reaching  the  observer 
at  azimuth  AZ,  and 

AZ 

is  the  instrument  azimuth  angle  relative  to  solar 
azimuth  (defined  as  zero  when  it  is  directly  away 
from  the  sun) . 

3.4.4  Equational  Form 

i of  Scattering  Function 

After  the  data  were  corrected  for  path  length  and  time,  the  azimuth  de- 
pendent radiance  data  were  normalized  by  the  data  point  nearest  to  pure 


* See  Section  2 of  this  volume,  where  normalization  of  this  measurement 
is  discussed. 
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backscatter  and  the  relative  azimuth  was  redefined  as  the  angle  from  back- 
scatter,  #.  The  model,  fitted  to  this  normalized  data  was: 

f(#)  = Aq  + Aj(l  + cos2#)  + cos2  Equation  3-13 

where  AQ,  A^,  A^,  A^  are  the  weighting  coefficients. 

The  terms  of  this  equation  were  arrived  at  by  examination  of  the  normalized 
polar  scatter  function  measurements.  It  involved  adding  the  classical 
Rayleigh  term  to  an  empirical  approximation  of  the  aggregate  scatter  profile 
of  atmospheric  particles  on  the  range  of  Mie  theory  sizes. 

The  data  for  each  of  the  three  wavelengths  was  analyzed  separately  for  the 
four  coefficients.  Data  points  within  7 degrees  of  the  point  180  degrees 
from  backscatter,  where  the  instrument  is  pointed  directly  toward  the  sun, 
were  eliminated  because  of  the  influence  of  direct  sunlight  on  the  instrument 
measurement . 

For  each  data  set,  the  four  coefficients  (AQ,  A^,  A2>  A^)  were  found  by  use 
of  nonlinear  regression  routines.  For  each  wavelength,  the  coefficients 
were  further  analyzed  as  functions  of  moisture  scale  height  and  the  Aj  co- 
efficient. The  wavelength  dependency  was  then  incorporated  by  finding  the 
coefficients  for  the  moisture  scale  height  and  the  AQ  3 terms  as  log- linear 
fits  to  the  450-nanometer,  650-nanometer,  and  850-nanometer  wavelength  mea- 
surements. The  final  equations  for  computing  the  A^,  A^,  A2>  and  A3  co- 
efficients of  Equation  3-13  were: 

Aq  = 0.765  - 1.412A1  Equation  3-14 

-2  1 

A:  = 0.947e  - (0.258  + 0.243  log1QX)  MSH  Equation  3-15 

A2  = 0.288  - 0.813A1  , and 


Equation  3-16 


8.083  ♦ 7.076  log1()A  - (10.37  ♦ 7.836  log1()A)  Aj 


Equation  3-17 


where  X is  the  wavelength  in  nanometers  and  MSH  is  the 
moisture  scale  height  in  centimeters. 


3.4.5  Volumetric  Normalization 

The  function  in  Equation  3-13  was  then  volumetrically  normalized  to  corre- 
spond with  the  definition  of  the  scattering  extinction  coefficient.  To 
accomplish  this,  it  was  necessary  to  find  the  scattering  volume  by  integrat- 
ing the  scattering  profile  over  all  solid  angles: 


Equation  3-18 


2a  a 

/«♦>  " ■ //  f(*)  sin*  d*  de 

n oo 

Where:  4 is  the  angle  from  backscatter  in  the  plane  of  the 

sun's  rays  and 

8 is  the  angle  of  rotation  about  the  sun's  ray. 


Substitution  of  Equation  3-13  into  Equation  3-18  results  in  the  following 
equation: 

n/2  a/2 

ft  (♦)  dfl  ■ 4a  Aq  f sin*  d*  ♦ 4a  A ^ J (1  + cos2*)  sin  4 d* 


a/3 


4a  A2  J cos2( sin*  d* 

0 

4a  Aj  sin*  d* 


Equation  3-19 
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Because  the  terms  of  Equation  3-19  could  not  all  be  solved  in  closed 
form,  numerical  integration  techniques  were  used.  The  resulting  inte- 
gral was  a function  of  the  coefficients  An  * . 


J f(4>)dn  * 4nA0  + 5.33uAj  ♦ .938ttA2  + .315irA3 


Equation  3-20 


Dividing  Equation  3-13  through  by  the  normalizing  integral  of  Equation  3-20 
results  in  polor  scatter  function**  that  was  coded  into  the  program  SCAT3, 
as  follows: 


/f(*)  = Bq  ♦ B1  (1  ♦ cos2*)  + B2cos2(^-)  ♦ B3e"3U'1t)‘ 


2J  3 


Equation  3-21 


Where:  B~  = -r- 
0 4tt 


B1  ‘ 4¥ 


1 / ^1  ^2 

-T*—[  1 + 1.333  — + .235  — + .079  ) 

4m  ^ A0  Aq  A0/ 

1 / An  ^3  \ 

-f—  ( — + 1.333  + .235  -r=—  + .079  -r^~  \ 
4,  ^ Aj  Aj  Aj  j 

4 if  (4;  * 1-333  4j‘  -23S  * -079  sr) 


B,  . ♦ 1.333  -1-  . .235  * .079  \ , 

3 «•  y ‘3  *3  ) 

Aq,  Aj,  A2,  and  A3  are  computed  from  Equations  3-14 

through  3-17. 


Equation  3-21  is  useful  because  it  always  integrates,  over  all  solid  angles, 
to  a unit  volume  regardless  of  the  estimates  of  AQ,  Aj,  A2>  and  Aj. 


* The  A.  coefficient  is  exactly  equal  to  the  w integral  for  pure  Rayleigh 
scatter.  3 

**  Definition  of  this  function  was  also  necessary  in  the  analysis  of  atmo- 
spheric path  spectral  radiance  as  discussed  above  in  para.  3.1. 
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3.4.6  Examples  of  Scatter  Function  Fit  to  Data 


Figure  3-9  shows  typical  examples  of  the  model  expressed  in  Equation  3-21, 
fitted  to  typical  examples  of  the  original  polar  scatter  data  corrected  for 
time  and  path  length. 

3.5  SCAT 3 MODEL  FITTED  TO  MEASURED  DATA 

Equation  3-12  was  used  to  generate  spectral  atmospheric  path  radiance  func- 
tions corresponding  to  specific  measurements  on  days  typical  of  the  collection 
program.  The  actual  measured  data  and  the  model  fit  are  shown  in  Appendix  C, 
Figures  C-l  through  C-29. 
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SECTION  4 

ANALYSIS  AND  MODELING  OF  HORIZONTAL  DAYLIGHT  IRRADIANCE 


4.1  APPROACH 

The  total  irradiance  incident  upon  a horizontal  surface  can  be  considered 
to  be  the  sum  of  three  distinct  irradiance  forms:  (1)  direct  solar  irradia- 
tion, (2)  indirect  atmospheric  downward  scatter,  and  (3)  indirect  reflected 
irradiance  from  the  earth's  surface  rescattered  downward  by  the  atmosphere. 

To  create  the  spectral  model  described  below,  each  individual  horizontal  day- 
light irradiance  measurement  made  by  the  radiometric  trailer  was  separated 


In  preparing  the  data  set  for  analysis,  only  those  day-sites  having 
estimates  of  the  spectral  transmittance  parameter  (see  Section  2 of  this 
volume)  were  used.  This  eliminated  most  of  the  days  having  obvious  cloud 
structure  and  permitted  the  analysis  to  be  devoted  to  modeling  of  uniform 
but  optically  varying  cloud- free  atmospheric  conditions.  Data  from  days 
having  cloud  cover  observable  from  the  ground  were  only  included  if  there 
was  no  indication  of  sun  obscuration. 

4.2  REGRESSION  ANALYSIS 

Expressed  mathematically,  the  component  sum  equivalent  to  total  horizontal 
daylight  irradiance  according  to  Robinson17  is: 

HDH(X)  = HDCX)  + bsHS(X)  + WX)  Equation  4-1 

Where:  H H(X)  is  the  total  measured  horizontal  daylight  spectral 
irradiance  in  watts/m2/5nm 

H (X)  is  the  direct  solar  spectral  irradiance  in  watts/ 
m2/5nm 

H„(X)  is  the  indirect  skylight  spectral  irradiance  in 
watts/m2/5nm 

H.  (X)  is  the  indirect  skylight  spectral  irradiance  derived 
from  reflected  energy  from  the  earth  (albedo  term) 
in  watts/m2/5nm,  and 

b and  bfl  are  the  proportionality  constants  derived  from  re- 
s gression  analysis  of  the  measured  data. 

4.2.1  The  Direct  Solar  Irradiance  Term 

The  first  term,  HD(X),  can  be  described  explicitly  in  terms  of  path  transmit- 
tance, the  cosine  of  the  solar  zenith  angle,  and  the  extraterrestrial  irradi- 
ance as  shown  below. 

VX)  = HSC(X)  * C0S6  IVX>  * ta(X)  ]M(6)  Equation  4-2 

17  N.  Robinson,  Solar  Radiation,  Solar  Physics  Laboratory,  Haifa,  Israel, 

1966,  pp  118-120. 
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Where:  HgC(X)  is  the  extraterrestrial  spectral  irradiance 
(solar  constant)  in  watts/m2/S  nm 

0 is  the  solar  zenith  angle 

t (X)  is  the  vertical  spectral  transmittance  due  to 
s scattering  phenomena 

ta(X)  is  the  vertical  spectral  transmittance  due  to  absorption 
phenomena,  and 

M(9)  is  the  air  mass  corresponding  to  zenith  angle  0. 

4.2.2  The  Indirect  Skylight  Irradiance  Term 

The  definitions  of  the  second  two  terms  of  Equation  4-1,  Hg(X)  and  HA(X) , are 
less  easy  to  derive  because  they  involve  the  integration  of  the  non-uniform 
downward  sky  hemisphere  irradiance,  and  in  the  case  of  H^Cx) , the  non-uniform 
surround  reflectance  of  the  terrain.  The  equational  forms  chosen  for  use  in 
the  model  were: 

H„(X)  = Hcr(X)  • cos0  [1  - t M(6)(X)]  • t (X)  Equation  4-3 

D oL  S cl 

and 

Ha(X)  = HgC(X)  • cos  [ta(X)  • tg(X)  ]MfX)R(X)  • Equation  4-4 

[1  - ts(X)]  • ta(X) 

where  R(X)  is  the  spectral  reflectance  of  the  general  surround  terrain.* 

Substituting  Equations  4-2,  4-3,  and  4-4  into  Equation  4-1,  dividing  through 
by  common  terms,  and  subtracting  1.0  from  both  sides  results  in: 


* See  Volume  2,  Section  5.1  of  this  report  for  the  derivation  of  these 
reflectance  data. 
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"dh(>) 

H (X)  • cose  Its(X)ta(X)] 


m(0)  ^ bs 


[i-ts(x)]m(e)ta(e) 


[ta(x)ts(x)]‘ 


Equation  4-5 


4 b R(X)  [1-t  (X)J  t (X) 


Examination  of  the  measured  spectral  irradiance  data  showed  that  the  terrain 
term  did  not  noticeably  affect  the  irradiance  data  except  under  extreme 
conditions  of  snow  cover.  By  disregarding  the  day-site  measurements  where 
snow  was  present,  the  bg  coefficient  was  found  by  linear  regression  tech- 
niques for  all  other  spectral  samrleo.  The  spectral  transmittance  values 
t (X)  and  t (X)  were  derived  from  ».ne  scattering  and  absorption  coefficients 
estimated  for  each  day  of  measurement*. 


Within  the  random  error  of  estimation,  the  daily  plots  of  bg  versus  solar 
altitude  could  be  estimated  from  linear  equations.  Figure  4-2  shows  two 
typical  examples  of  the  bg  coefficient  versus  solar  altitude.  The  slope  and 
intercept  of  the  daily  plots  of  bg  were  regressed  against  moisture  scale 
height  and  found  to  correlate.  Figures  4-3  and  4-4  show  the  slope  and  inter- 
cept terms  versus  moisture  scale  height  and  fitted  with  the  regression  equa- 
tions. The  complete  equation  for  the  indirect  skylight  irradiance  coefficient 
bs  in  terms  of  solar  altitude  and  moisture  scale  height  is: 


bg  = .387  - . 0131MSH  - (.00392  - .00221MSH)  SA 


Where:  MSH  is  moisture  scale  height  (cm),  and 
SA  is  solar  altitude  (degrees) 


Equation  4-6 


See  Section  2,  Table  2-1,  of  this  volume  for  the  list  of  coefficients. 


Bs  COEFFICIENT 


DATE:  6/11/76 
MOISTURE  SCAI£  HEIGHT 


SOLAR  ALTITUDE 


DATE:  7/21/76 

MOISTURE  SCALE  HEIGHT  - 3.03  cm 


MODEL  

ACTUAL  DATA  X 


SOLAR  ALTITUDE 


Figure  4-2.  Two  Typical  Plots  Showing  the  Bs  Coefficient 
as  a Function  of  Solar  Altitude 


COEFFICIENT  OF  THE  LINEAR 
SOLAR  ALTITUDE  TERM* 


where  the  equation  limits  are:  0 < MSH  < 3.5  cm 

-5°  < SA  < 90° 


4.2.3  Terrain  Reflectance  Term 


The  day-site  data  collections  in  which  there  was  extensive  snow  cover  were 
used  to  estimate  the  coefficient  of  the  terrain  reflectance  term.  By  evalu- 
ating the  bg  coefficient  for  the  appropriate  solar  altitude  and  moisture 
scale  height,  the  skylight  term  could  be  effectively  subtracted  from  the 
measurement  in  the  same  manor  as  the  direct  solar  irradiance  term  and  resulted 
in  the  following  equation: 

Hdh(A)  [1-t  (A)m(e)]  t (A) 

;57eT  -l  -b  (SA.MSH)  5 * 

Hsc(A)  • cos0[ts(A)ta(A)mte-1  [ts(A)  ta(A)]mCe) 

Equation  4-7 

= baR(A)  [1— 1$ (A) ] ta(A) 


By  following  the  same  regression  procedure  used  for  b , b coefficients  for 
each  solar  altitude  and  each  day-site  could  be  estimated.  A solar  altitude 
plot  of  the  average  b coefficient  for  all  snow  sites  is  shown  in  Figure  4-5 

a 

The  upward  curvature  of  ba  at  low  solar  altitude  was  considered  to  be  more 
an  artifact  of  the  analysis  than  a physical  condition  and  was  therefore  not 
modeled.  No  attempt  to  correlate  the  bfl  coefficient  with  moisture  scale 
height  was,  or  could  be,  made.  The  resulting  equation  for  b was  therefore: 

3 


b = 3.908  - 0.0739  SA 

a 


Equation  4-8 


The  value  of  bfl  is  limited  to  a minimum  value  of  0.5  corresponding  approx- 
imately to  the  highest  solar  altitude  observed  with  snow  surround. 


4.3  MODEL  EQUATIONS  IN  SCAT 3 

The  mathematical  model  for  horizontal  daylight  irradiance  discussed  in  this 
section  is  a mixture  of  physical  and  empirical  equations.  It  has  the  dis- 
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b pi  f 

SOLAR  ALTITUDE 


Figure  4-5.  Average  Terrain  Term  Coefficient 
as  a Function  of  Solar  Altitude 


tinct  advantage  of  satisfying  both  the  requirements  of  atmospheric  physics 
and  observed  spectral  data.  In  the  final  form  of  the  SCAT3  equation,  the  spec- 
tral reconstruction  of  horizontal  daylight  is  accomplished  as  follows: 


hdhU)  = hsc(x)  cose  [ts(x)ta(x)]m(e)  . kD 


♦ kg  • bs(SA,MSH) 

* ba(5A)  • Hsc(x) 


HSC(X)  cose 


ta(x) 


cose  • [ts(x)ta(x)]m(e) 


Equation  4-9 


• R(X)  [l-ts(x)]  ta(x) 


Where:  Hg„(X)  is  the  mean  extraterrestrial  spectral  solar 
irradiance  corrected  for  seasonal  variation* 
(watts/m2/5nm) , 

0 is  the  solar  zenith  angle  between  0°  and  87°, 

t (A),t  (A)  are  the  spectral  vertical  atmospheric  scatter- 
ing and  absorption  transmittances, 

M(e)  is  the  relative  air  mass  corresponding  to  solar 
zenith  angle  6, 

RCA)  is  the  spectral  reflectance  of  the  general  terrain, 

bs(SA,MSH)  is  the  skylight  irradiance  proportionality  constant 

evaluated  for  a specific  solar  altitude  and  mois- 
ture scale  height  (as  calculated  in  Equation  4-6), 

b (SA)  is  terrain/skylight  irradiance  proportionality 

constant  evaluated  for  a specific  solar  altitude 
(Equation  4-8),  and 


* See  Appendix  D for  a description  of  this  correction  which  adjusts  the 
solar  constant  for  the  orbital  eccentricity  of  the  earth  about  the  sun. 
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. ) 1,  SA  > 0° 

kD  * j 0.  SA  < 0° 

. ) 0,  SA  > 3° 

s “ ) p(SA)  SA  < 3°* 

SA  = Apparent  solar  altitude** 


4.4  VERY  LOW  SOLAR  ALTITUDE  MODELING 

Only  a few  daylight  irradiance  measurements  were  available  below  three 
degrees  solar  altitude.  To  continue  the  model's  applicability  down  to 
minus  five  degrees  solar  altitude,  it  was  necessary  to  adopt  a relative 
scaling  polynomial  for  daylight  that  was  derived  from  separate  data. 

The  polynomial  itself  does  not  involve  spectral  shaping;  in  addition,  it 
contains  no  parameters  other  than  solar  altitude  and  is  used  strictly  as 
a relative  magnitude  adjustment  for  the  skylight  irradiance  term,  Hg(X). 

The  adjustment  polynomial  p(SA)  had  been  evaluated  at  three  degrees  solar 
altitude,  and  this  value  was  divided  into  the  coefficients  to  normalize 
the  polynomial  to  the  lowest  valid  calculated  daylight  spectral  irradiance. 
Equations  4-1  through  4-4  for  the  model  are  evaluated  for  the  appropriate 
condition  of  moisture  scale  height  and  terrain  reflectance  at  three  degrees 
solar  altitude;  they  are  then  multiplied  wavelength-by-wavelength  by  the 
normalized  adjustment  polynomial  evaluated  at  the  solar  altitude  of  inter- 
est. The  polynomial  and  scaling  approach  are  explained  in  Appendix  F. 

The  normalized  factor  used  for  horizontal  daylight  spectral  irradiance  is 
0.037440. 

Because  refraction  of  the  sun's  rays  in  the  atmosphere  becomes  substantial 
below  five  degrees  solar  altitude,  a correction  must  be  applied  to  solar 
altitude  as  was  used  in  Equation  4-10.  The  corrected  solar  altitude  is  re- 
ferred to  as  the  "apparent"  solar  altitude  and  is  described  mathematically 
in  Appendix  G. 


* See  para.  4.4  and  Appendix  F for  a description  of  the  polynomial 
P(SA). 

**  See  Appendix  G for  the  correction  applied  to  solar  altitude  below  5° 
to  compensate  for  refraction  of  the  sun's  rays  in  the  atmosphere. 
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4.5  SPECTRAL  RECONSTRUCTION  OF  MEASURED  DATA 


Equation  4-10  for  the  model  was  used  to  reconstruct  a series  of  actual 
spectral  daylight  irradiance  measurements  that  typify  the  overall  collec- 
tion program.  These  data  and  the  model  fits  are  shown  in  Appendix  C. 





SECTION  5 

ANALYSIS  AND  MODELING  OF  HORIZONTAL  SKYLIGHT  IRRADIANCE 

5.1  INTRODUCTION 

For  the  purposes  of  this  study,  horizontal  skylight  irradiance  is  defined  as 
the  total  irradiance  on  a horizontal  surface  with  the  direct  rays  of  the  sun 
blocked.  The  horizontal  surface  in  this  case  can  be  considered  to  be  irradi- 
ated by  light  from  the  sun  that  is  both  (1)  scattered  downward  by  the  atmo- 
sphere on  the  horizontal  surface  and  (2)  reflected  from  the  surrounding  terrain 
and  scattered  by  the  atmosphere  downward  to  the  horizontal  surface. 

To  create  the  model  for  horizontal  skylight  irradiance,  the  following  were 
used  in  the  analysis:  irradiometer  data  recorded  by  the  mobile  radiometric 
trailer,  estimates  of  the  vertical  transmittance  based  on  the  analysis  in  Sec- 
tion 2,  and  meteorological  data  collected  concurrent  with  the  radiometric  data. 
The  data  had  a solar  altitude  range  of  0.73  to  80.95  degrees  and  a moisture 
scale  height  range  of  0.23  to  4.2  c _*nt imeters . Figure  5-1  illustrates  the 
geometry  of  the  model  used  in  this  analysis.  Each  of  the  additive  terms  has 
a coefficient  associated  with  it  that  was  estimated  by  linear  regression  of 
the  model  equational  terms  to  the  measured  skylight  irradiance. 


Figure  5-1.  Geometry  Model  for  Horizontal  Skylight  Irradiance 
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The  skylight  irradiance  model  is  illustrated  in  Figure  5-1  and  shown  mathe- 
matically below. 

HSH(X)  " bSHS(X)  + WX)  Equation  5-1 


Where:  H (X) 

is  the  total  horizontal  skylight  spectral 
irradiance,  in  watts/m2/Snm 

HS(A) 

is  the  spectral  irradiance  from  the  atmospheric  scatter 
of  direct  solar  rays,  in  watts/m2/5nm 

VX) 

is  the  spectral  irradiance  from  atmospheric  scatter 
of  reflected  daylight  from  the  surround  terrain 
reflectance,  in  watts/m2/5nm. 

bS,bA 

are  the  proportionality  constants  determined  from  re- 
gression analysis. 

Derivation  of  the  above  terms  and  their  coefficients  is  explained  below. 

5.2  ANALYSIS  OF  COMPONENTS 

As  was  the  case  for  other  forms  of  ground- level  irradiance,  the  equation 
terms  were  taken  sequentially  in  the  regression  analysis. 

5.2.1  Direct  Atmospheric  Scatter  Term 

The  first  component  of  skylight  irradiance  is  the  direct  solar  irradiance 
scattered  downward.  This  term  represents  a fraction  of  the  total  irradiance 
outside  the  atmosphere  that  is  not  transmitted  but  scattered  by  the  atmo- 
sphere, and  is  represented  by  the  equation: 

Hs(x)  " Hsc(x)cos0  * [*  - TS(e)(^]‘  Vx)  Equation  5-2 

Where:  HS„(A)  is  the  extraterrestrial  spectral  irradiance,  in 
watts/m2/5nm, 

T_(A)  is  the  vertical  spectral  transmittance  due  to  scatter- 
ing. 

T.  (A)  is  the  vertical  spectral  transmittance  due  to  absorp- 
tion, and 

m(e)  is  the  air  mass  at  solar  zenith  angle  6. 
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It  was  assumed  that  the  terrain  reflectance  term  of  Equation  5-1,  HA(A), 
was  negligible  for  all  day-sites  except  those  with  snow  surround,  and  all 
non-snow  day-site  data  was  regressed  according  to  the  following  equation: 


HSH^ 


= br 


HgC(A)cose 


[‘  - vtft6) 


ta(») 


Equation  5-3 


The  regression-determined  value  for  bg  was  then  studied  as  a function  of 
solar  altitude  and  moisture  scale  height.  Each  day-site  b value  set  was 
fitted  with  a linear  regression,  and  the  slope  and  intercept  were  determined. 
Figures  5-2  and  5-3  are  two  typical  plots  of  the  b„  coefficient  versus  solar 
altitude. 

No  relationship  with  moisture  scale  height  could  be  found  for  the  intercept 
of  the  linear  equations;  however,  as  Figure  5-4  illustrates,  the  slope  term 
correlated  with  moisture  scale  height  to  a slight  degree.  The  final  equation 
for  the  b,,  coefficient  was  therefore: 

bg(SA,MSH)  = 0.348  + (-0.000326  + 0.001318  MSH)  SA  Equation  5-4 

Where:  MSH  is  the  moisture  scale  height,  between 
0 and  3.5cm,  and 

SA  is  solar  altitude  (SA  = 90°  - G) . 
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Figure  5-2.  Typical  Example  of  the  bs  Coefficient  Relationship  to 
Solar  Altitude  for  Low  Moisture  Scale  Height 


Figure  5-3.  Typical  Example  of  the  bs  Coefficient  Relationship  to 
Solar  Altitude  for  High  Moisture  Scale  Height 
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5.2.2  Terrain  Reflectance  Term 


The  coefficient  for  the  terrain  term,  b^,  was  investigated  using  day- 
site  measurements  acquired  while  snow  was  present  on  the  ground  around  the 
measurement  site.  An  equational  form  identical  with  that  for  horizontal 
daylight  irradiance  was  chosen  for  use  and  regressed  against  the  measured 
data  as  follows: 

[l-ts(A)] 

bs  ^ [l-tsm(e)(A)] 

[t  (A).  t fA)]m(0)-  R(A)  Equation  5-5 

a S 

where  R(A)  is  the  spectral  reflectance  of  the  surrounding  terrain. 


hsh(x) 

hsc(>0  cose  [i-ts"te)(x)]  ta(J) 


The  terrain  coefficient,  bA,  was  found  to  operate  in  a manner  similar  to 
that  for  horizontal  daylight  irradiance  (see  para.  4.2  of  this  volume). 

The  horizontal  daylight  functions  for  bA  below  were  therefore  used  for  the 
skylight  irradiance  model. 

Equation  5-6 

bA(SA)  = 3.908  - 0.0739  • SA 
where  SA  is  the  solar  altitude  (SA  = 90°  - 6) . 


5.3  VERY  LOW  SOLAR  ALTITUDE  MODELING 

The  model  described  in  Section  5 was  considered  to  be  valid  to  a solar  alti- 
tude as  low  as  3 degrees.  Extension  of  the  model  to  -5  degrees  solar  altitude 
was  accomplished  using  the  scaling  function  described  for  the  horizontal  day- 
light irradiance  model,  in  Section  4,  and  Appendix  F of  this  volume. 

The  scaling  factor  is  applied  to  the  skylight  term  and  operates  only  below 
3 degrees  solar  altitude.  The  factor  is  not  a function  of  wavelength  and 
atmospheric  condition. 


5.4  MODEL  EQUATIONS  IN  SCAT3 


The  following  equation  is  used  in  the  program  to  model  horizontal  skylight 
irradiance: 


HSH(X)  = ks  * bS  (SA’MSH)  * HSC(A)  cos6  [l-^smC0) (A)]  • tfl(X)  Equation  5-7 
♦ bAHsc(X)  cose  [ta(X)ts(X)]m(0).  R(X) 


[l-ts(X)]  • ta(X) 

Where:  Hgc(x)  is  the  mean  extraterrestrial  spectral  solar 
irradiance  corrected  for  seasonal  variation* 
(watts/m2/5nm) . 

ts(*),ta(X)  are  the  vertical  spectral  transmittances  due 
to  scattering  and  absorption. 


R(X)  is  spectral  reflectance  of  terrain. 

M(e)  is  the  relative  air  mass  corresponding  to  the 

solar  zenith  angle  0. 

0 is  the  solar  zenith  angle,  between  08  and  87°. 


b (SA,MSH)  is  the  skylight  irradiance  proportionality  con- 
stant evaluated  for  a specific  solar  altitude 
and  moisture  scale  height  (as  computed  in  Equa- 
tion 5-4) . 


bA(SA) 

SA 


is  the  terrain/skylight  irradiance  proportionality 
constant  evaluated  for  a specific  solar  altitude 
(as  computed  in  Equation  5-6). 


is  the  apparent  solar  altitude,**  and 


1,  SA  * 3° 

p(SA) , SA  < 38*** . 


* See  Appendix  I 

**  See  Appendix  G 

See  Appendix  F 


*** 


5.5  SPECTRAL  RECONSTRUCTION  OF  MEASURED  DATA 


Equations  5-4,  5-6,  and  5-7  were  used  to  re-create  spectral  skylight  irradi- 
ance  data  corresponding  to  a typical  set  of  measured  data.  Graphs  of  these 
data  are  shown  in  Appendix  C. 


SECTION  6 

ANALYSIS  AND  MODELING  OF  VERTICAL  DAYLIGHT  IRRADIANCE 
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6.1  INTRODUCTION 

The  vertical  surface  irradiance  measurements  made  by  the  mobile  radio- 
metric  trailer  and  described  in  Volume  2 of  this  report  were  all  made  with 
the  sampling  aperture  axis  co-planar  with  the  sun's  azimuth.  These  measure- 
ments were  of  two  types:  (1)  looking  directly  toward  the  sun,  and  (2)  look- 
ing directly  away  from  the  sun.  The  analysis  of  data  from  the  first  type, 
vertical  daylight  irradiance,  is  described  in  this  section.  The  second 
type  is  described  in  Section  7 of  this  volume. 

Although  the  terms  are  similar,  the  model  for  vertical  front  daylight  irradi- 
ance differs  in  several  respects  from  that  for  horizontal  daylight  as  shown 
below. 


HDV(A>  = 


H'CA) 


+ b' 
s 


Hg(A) 


H^(A) 


Equation  6-1 


Where:  HDV(A)  is  the  spectral  vertical  daylight  irradiance 
(front-lighted),  in  watts/m2/5nm, 


H£(A) 


is  the  direct  solar  spectral  irradiance,  in 
watts/m2/5nm. 


H'(X)  is  the  indirect  skylight  irradiance,  in 
watts/m2/5nm) . 


!1'(X)  is  the  indirect  reflected  irradiance  from  the  nearby 
a terrain,  in  watts/m2/5nm,  and 

b',b'  are  the  proportionality  constants  determined  from 
s a regression  analysis. 
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Figure  6-1  shows  the  geometry  that  is  assumed  in  the  basic  model  of 
Equation  6-1. 


Figure  6-1.  Geometry  Model  Used  for  Vertical  Daylight  Irradiance 


6.2  DIRECT  SOLAR  IRRADIANCE  COMPONENT 

The  direct  sunlight  component  can  be  expressed  as  a product  of  four  items: 
the  solar  constant  H (X)  times  the  scattering  transmittance  t (X)m^ 
times  the  absorption  transmittance  t (X)  ^ ' times  a sine  of  the  solar  zenith 
angle  which  is  in  conformity  with  Lambert's  law  of  illumination.  Thus, 
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the  irradiance  of  the  direct  component  can  be  written  as: 

Hp(X)  = Hg^(X)  sin©  • [ts(X)  • ta(X)]m^  Equation  6-2 

Where:  Hg^X)  is  the  spectral  extraterrestrial  solar  irradi- 
ance in  watts/m2/ 5nm. 

ts(X),ta(X)  are  the  scattering  and  absorption  spectral 
transmittances . 

m(6)  is  the  air  mass  corresponding  tc  che  solar  zenith 
angle  6,  and 

9 is  the  solar  zenith  angle  (equal  to  90°  solar  altitude) . 


6.3  INDIRECT  SKYLIGHT  COMPONENT 

An  equation  that  approximates  the  skylight  term  can  be  formed  by  considering 
the  atmosphere  as  a plane  surface  that  is  irradiated  by  the  sun.  The  light 
that  passes  through  the  atmosphere  is  partially  scattered  and  partially 
absorbed.  The  scattered  portion  is(l-Tg'n^6^)  and  the  portion  not  absorbed 
is  Ta.  Thus,  the  component  of  the  total  skylight  reaching  a vertical  surface 
can  be  approximated  by: 

H'(X)  = Hc  • cos0  • [1-t _(X)m^]  • t (X)  Equation  6-3 

O OL.  S 3 

6.4  INDIRECT  SURROUND  REFLECTANCE  COMPONENT 

Since  any  plane  orientation  other  than  horizontal  receives  some  irradiance 
from  the  terrain  within  the  field  of  view,  a third  term  is  necessary  that 
involves  the  reflectance  of  the  instrument  surround.  The  mathematical  approx- 
imation for  this  term  is: 

H'(X)  - Hcr(X)  • cos0  • [t  (X)  • t (X)]m(X)  • R(X)  Equation  6-4 

3 oL*  S 3 
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where  R(X)  is  the  spectral  reflectance  of  the  terrain  surrounding  the 
instrument.* 

6.5  REGRESSION  ANALYSIS  FOR  PROPORTIONALITY  CONSTANTS 

The  three  terms  described  in  Equations  6-2,  6-3,  and  6-4  were  substituted 
into  Equation  6-1  and  algebraicly  manipulated  as  follows: 


HPV^ 

HSC(A)  cose  [ta(X)  • ts(X)]m(6) 


tane  = b' 

S 


[:-ts(x)m(0)]  ta(x) 

[ta(A)  • ts(A)]mC0) 


Equation  6-5 


+ b;  R(A)[l-ts(A)]  • ta(X) 


Observations  that  had  been  made  from  the  measured  data  indicated  that,  as 
in  the  case  of  horizontal  daylight  irradiance,  the  indirect  skylight  com- 
ponent coefficient  b^  could  be  determined  from  regression  analysis  of  the 
data  from  all  days  not  having  snow  surround  so  that  the  terrain  reflectance 
term  could  be  assumed  to  be  zero.  The  b^  coefficient  was  then  studied  as  a 
function  of  solar  altitude  and  moisture  scale  height.  The  coefficient  did 
not  show  strong  correlation  with  either  solar  altitude  or  moisture  scale 
height.  However,  a linear  model  was  regressed  and  fitted  because  of  a weak 
dependence  upon  these  parameters. 

The  assumption  that  the  terrain  reflectance  term  is  small  compared  to  the 
skylight  contribution  for  low  and  medium  surround  reflectances  was  verified 
by  experiment.  Vertical  daylight  irradiance  measurements  were  made  in  the 
normal  manner  that  were  followed  immediately  by  similar  measurements  made  with 
a black  cloth  positioned  so  that  it  eliminated  direct  reflections  from  the 


* Section  5.1,  Volume  2,  for  a description  of  these  data. 


terrain  into  the  irradiometer.  These  measurements,  when  compared  with 
the  direct  component  calculated  from  the  measured  transmittance,  showed 
the  skylight  contribution  to  be  considerably  greater  than  that  from  the 
terrain  reflectance. 

6.5.1  Skylight  Coefficient 

The  bg  coefficient  that  corresponds  to  the  proportion  of  skylight  irradi- 
ance,  is  shown  in  Figure  6-2  as  a function  of  the  moisture  scale  height. 

The  plotted  bg  coefficient,  which  is  the  average  coefficient  determined  for 
all  solar  altitudes  of  a given  day-site,  shows  a tendency  to  increase  with 
moisture  scale  height.  A slight  solar  altitude  dependence  of  bg  can  be  seen 
in  Figure  6-3,  which  is  a graph  of  all  the  values  of  bg  for  each  day-site 
plotted  versus  their  respective  solar  altitudes.  The  regression  equation 
shown  in  Figure  6-2  for  bg  versus  moisture  scale  height  was  determined  first. 

A multiplying  factor  was  then  determined  that  was  based  on  the  solar  altitude 
versus  b regression,  and  it  was  used  to  modify  the  b moisture  scale  height 
relationship. 

The  multiplying  factor  is  the  value  of  bg  at  a given  solar  altitude  as  deter- 
mined by  the  regression  equation  shown  in  Figure  6-3  divided  by  the  regression 
value  of  bg  at  28  degrees  solar  altitude  (the  median  solar  altitude  of  the 
data  forming  the  relationship  of  Figure  6-2).  Thus,  the  multiplying  factor  is 
normalized  so  that  it  does  not  change  bg  at  28  degrees  solar  altitude,  but 
makes  a change  at  all  other  solar  altitudes  equal  to  the  fractional  change  of 
bg  in  the  solar-altitude-dependent  regression  of  Figure  6-3.  The  equation  for 
bg  in  terms  of  both  the  solar  altitude  factor  and  the  moisture  scale  height 
becomes: 

bg  = (0.3355  + 0.0428  x MSH)  x (0.848  + 0.00544  x SA)  Equation  6-6 

where  MSH  is  the  moisture  scale  height  between  0 and  4.2cm,  and 
SA  is  the  solar  altitude,  in  degrees. 


MOISTURE  SCALE  HEIGHT  (cm) 


Figure  6-2.  Vertical  Front-Lit  Skylight  Coefficient  as  a Function 
of  Moisture  Scale  Height 
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6.5.2  Terrain  Reflectance  Coefficient 

A scaling  coefficient,  bA>  that  was  based  on  acquired  data  was  calculated 
for  the  term  expressed  in  Equation  6-4.  Because  only  the  days  with  snow 
surround  were  considered  to  have  sufficient  reflectance  to  exceed  the  noise 
level  of  the  measurements,  only  these  data  were  considered  in  forming  the 
bA  coefficient.  In  Figure  6-4,  it  can  be  seen  that  bA  is  a function  of 
solar  altitude.  At  about  10  degrees  solar  altitude,  specular  reflections 
caused  a large  increase  in  the  bA  coefficient,  making  it  necessary  to  model 
the  bA  coefficient  in  two  parts  (below  10°  SA  and  above  10°  SA) . Both  solar 
altitude  ranges  were  fitted  separately  with  straight  lines  determined  by  least 
squares  regressions.  The  two  line  segments  were  then  modified  so  that  they 
would  be  without  discontinuity  at  10  degrees  solar  altitude.  The  final 
equation  for  bA  was  then: 


2.50  - 0.1945  SA,  SA  < 10° 
0.5927  - 0.003769  SA,  SA  > 10° 


Equation  6-7 


The  terrain  term  above  is  allowed  to  go  to  zero  at  zero  degrees  solar  altitude 
in  accordance  with  the  cosine  term  in  Equation  6-4.  It  is  set  to  zero  for 
solar  altitudes  less  than  zero  degrees. 


6.6  MODIFICATION  OF  THE  DIRECT  AND  SKYLIGHT  COMPONENTS  FOR  IMPROVED 
FITTING  TO  EMPIRICAL  DATA 


The  sum  of  the  direct,  skylight,  and  terrain  terms  as  calculated  in  the 
model  slightly  underestimated  the  vertical  daylight  irradiance  compared  to 
a cross-section  of  actual  measured  data.  To  make  them  fit  the  model  better, 
the  first  two  terms  (direct  and  skylight)  were  increased  by  a factor  calcu- 
lated to  correct  for  the  average  offset.  The  necessary  multiplying  factor 
was  found  to  be  1.14.  The  factor  was  not  applied  to  the  albedo  term,  how- 
ever, since  the  offset  was  needed  even  in  cases  where  the  terrain  contribu- 
tion could  be  considered  negligible. 
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6.7  VERY  LOW  SOLAR  ALTITUDE  MODELING 

Although  the  model's  equations  were  shown  to  be  valid  for  solar  altitudes 
down  to  and  including  6 degrees,  below  6 degrees,  there  was  insufficient 
data  to  verify  the  model.  Equation  6-2  continues  to  be  used  as  the  model 
for  the  direct  component  down  to  zero  degrees,  below  which  it  is  set  to 
zero.  However,  Equation  6-3  for  the  skylight  component  is  not  suitable 
down  to  zero  degrees  because  the  term  goes  to  zero  even  though  skylight  is 
physically  present  below  zero  degrees.  Modeling  of  skylight  below  6 de- 
grees is  accomplished  by  applying  the  adjustment  polynomial  P(SA)  given  in 
Appendix  F.  The  model-predicted  skylight  values  at  6 degrees  solar  alti- 
tude then  define  the  shape  of  the  function  with  the  skylight  irradiance 
decreasing  in  proportion  to  the  decrease  in  the  irradiance  scaling  polynomial 
p(SA) . 

6.8  EQUATIONS  USED  IN  SCAT3 

The  derivation  of  each  of  the  terms  making  up  the  vertical  daylight  irradi- 
ance were  discussed  above.  In  its  complete  form,  the  daylight  irradiance  is 
represented  by: 

HDV(A)  = kD  * HSC(A)  Sin6  * fts(A)  ta(A)1"(6) 

Equation  6-8 

+ 1.14  K*  • b'  (SA,MSH)  • Hcr  cos0  [1-t  rx)n,(6)  ] • t(X) 

S S 5 “ 

♦ b;  (sa)  • hsc  cose  [ts(x)  ta(x)]"(e)  • R(X) 

Where:  H (X)  is  the  spectral  extraterrestrial  solar  irradiance, 
in  watts/m2/5nm,  corrected  for  seasonal  variation*. 

6 is  the  solar  zenith  angle  (equal  to  90°  solar  alti- 
tude) . 

t (X),ta(X)  are  the  scattering  and  absorption  spectral  trans- 
mittance. 

"*  See  Appendix  I for  a description  of  this  correction. 
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* 


m(A)  is  the  air  mass  corresponding  to  the  solar  zenith  angle  6. 


R(9)  is  the  spectral  reflectance  of  the  terrain  surrounding 
the  instrument. 

b$ (SA.MSH)  is  the  skylight  irradiance  proportionality  constant  evalu- 
ated for  a specific  solar  altitude  and  moisture  scale  height 
computed  in  Equation  6-6. 

b^(SA)  is  the  terrain/skylight  irradiance  proportionality 

constant  evaluated  for  a specific  solar  altitude  com- 
puted in  Equation  6-7. 

SA  is  the  apparent  solar  altitude* 

K,  ) 1.  SA  > 0- 

( 0,  SA  < 0° 

K,  jl,  SA>6° 

s I P(SA),  SA  < 6°** 


* See  Appendix  G for  the  correction  applied  to  solar  altitude  below  5* 
to  compensate  for  refraction  of  the  sun's  rays  in  the  atmosphere. 

**  See  Section  4.4  and  Appendix  F for  a description  of  the  polynomial 
P(SA). 


SECTION  7 

VERTICAL  SKYLIGHT  IRRADIANCE 


7.1  INTRODUCTION 

Vertical  skylight  irradiance  is  defined  for  this  study  as  the  indirect 
irradiance  from  atmospheric  scatter  and  terrain  reflectance  on  a vertical 
plane  oriented  180  degrees  away  from  the  solar  azimuth.  This  measurement 
is  the  sum  of  contributions  from  the  skylight  half-hemisphere  and  the 
horizontal  surfaces  surrounding  the  instrument, that  define  the  terrain 
half-hemisphere.  Special  experiments  were  conducted  to  remove  the  effect 
of  surround  so  that  the  skylight  term  could  be  studied  independently.  The 
surround  term  was  later  re-introduced  into  the  final  equations.  Figure  7-1 
shows  the  model  assumed  in  the  analysis. 


Figure  7-1.  Geometry  Model  Used  for  Vertical  Skylight  Irradiance 
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7.2  ELIMINATION  OF  NEAR-FIELD  SURROUND  EFFECTS 


The  vertical  backlit  measurements  were  found  to  contain  a significant 
amount  of  trailer  reflection  irradiance  along  with  the  skylight  irradiance. 

The  amount  of  reflection  is  substantially  influenced  by  the  azimuth  of 
the  instrument  because  a variety  of  reflecting  surfaces  are  present  near 
the  irradiometers  on  the  trailer  roof.  As  an  irradiometer  changes  azimuth, 
it  sees  various  portions  of  the  ground,  the  trailer  top,  and  the  other 
instruments.  Furthermore,  the  white  styrofoam  hatch  covers  that  are  also 
visible  to  the  irradiometers  were  surfaces  of  higher  reflectance  than  the 
rest  of  the  gray  trailer  roof. 

In  order  to  back  the  unwanted  trailer  reflections  out  from  this  data,  an 
experiment  was  conducted  to  measure  vertical  backlit  skylight  at  various 
azimuths,  both  under  normal  operating  conditions  and  also  with  trailer  re- 
flections minimized.  The  trailer  reflections  were  eliminated  by  covering 
the  trailer  top  with  a large  black  canvas  tarp.  Measurements  were  also 
made  with  the  tarp  extended  from  the  trailer  so  as  to  block  almost  all 
direct  reflections  from  the  terrain  surrounding  the  trailer.  These  measure- 
ments indicated  that  ground  reflectances  were  insignificant  and  that  they 
could  generally  be  disregarded. 

Because  the  experiment  was  conducted  in  the  morning  of  two  consecutive  days, 
readings  were  made  with  Irradiometer  No.  1,  which  was  used  for  all  morning 
vertical  backlit  measurements  during  the  main  collection  program.  No  read- 
ings were  necessary  with  Irradiometer  No.  2,  which  had  been  used  for  after- 
noon vertical  backlit  measurements,  because  the  correction  for  its  measurements 
could  be  determined  from  symmetry  considerations. 

A plot  of  the  experimental  data,  both  with  and  without  trailer  reflections, 
is  shown  in  Figure  7-2  as  a function  of  azimuth  (550  nm  wavelength).  This 
plot  shows  that  an  increase  was  necessary  to  correct  the  data  as  computed 
from  the  ratio  of  the  readings  as  a function  of  azimuth.  It  was  also  observed 
that  the  correction  was  not  strongly  influenced  by  wavelength.  The  correction 
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255  260  265  270  275  281  287  293  300  309  318  329 

INSTRUMENT  AZIMUTH  RELATIVE  TO  NORTH  (IN  DEGREES) 

I . 

Figure  7-2.  Comparison  of  Vertical  Skylight  Irradiance  Data  With 
and  Without  Trailer  Reflections  (Values  at  550nm) 


ratio  was  computed  using  data  from  all  the  wavelengths  sampled  in  the 
experiment;  that  is,  550,  650,  750,  850,  950,  1050,  and  1150  nanometers. 

A third-order  polynomial  was  fitted  to  the  correction  ratios  as  follows: 

H'  (X)  = [3.52  - 0.307  A2  ♦ 0.000886  AZ2  - 0.000000823  AZ3] 

B Equation  7-1 

HSVB^ 

The  correction  ratios  for  Irradiometer  No.  2 were  made  by  substituting 
360°-AZ  for  AZ  in  Equation  7-1. 

7.3  REGRESSION  ANALYSIS 

Once  the  vertical  skylight  data  had  been  corrected  by  removing  the  effects 
of  the  trailer  roof,  a mathematical  model  was  fitted  for  regression  analysis. 

A surround  reflectance  term  was  later  added  by  assuming  that  the  surround 
contribution  would  be  the  high  solar  altitude  model  that  had  been  chosen 
for  vertical  daylight  (see  Section  6 of  this  volume). 

7.3.1  Skylight  Component 

The  skylight  can  be  approximately  modeled  by  considering  the  atmosphere  as 
a plane  surface,  the  illumination  of  which  is  equal  to  the  solar  constant 
times  the  cosine  of  the  zenith  angle  as  illustrated  in  Figure  7-1.  Taking 
absorption  and  scattering  into  account,  the  vertical  backlit  irradiance  can 
be  approximated  by  the  following  equation: 

IlVSK(X)  - b$  • Hsc(X)  cose  jl  - ts  V*)  Equation  7-2 

Where:  HVSK(>)  is  the  spectral  skylight  irradiance  component 

v K on  the  vertical  surface,  with  surround  reflec- 

tance effects  removed,  in  watts/m2/5nm. 

H (X)  is  the  mean  spectral  extraterrestrial  irradiance, 

sc  in  watts/m2/5nm. 


6 

is 

taU) 

is 

(A) 

is 

eff  scattering,  and 

bg  is  the  regression-determined  proportionality  constant. 

Unlike  the  other  regressions  discussed  above, the  model  had  to  be  modified 
to  better  approximate  the  spectral  shape  of  the  data.  The  correct  description 
for  t_  (X)  was  determined  by  trial  and  error,  first  by  trying  t (X)  * 


t$(A) 


eff 

m(6) 


eff 


This  equation  produced  skylight  values  that  were  too  high  in  the 
blue  end  of  the  spectrum  and  too  low  in  the  IR  and  of  the  spectrum.  A tailor- 
ing of  the  shape  of  this  function  was  achieved  by  letting  t (X)  equal  a 


fraction  of  t s ( X) 


m(0) 


as  in  the  following  equation: 
m(0) 


’eff 


*eff 


(X)  = 0.3  t (X) 


Equation  7-3 


Once  the  shape  had  been  defined,  the  scaling  coefficient  b could  be  com- 
puted to  make  the  skylight  model  fit  the  data.  This  scaling  coefficient,  as 
determined  by  the  regression  of  Equation  7-2,  is  a function  of  solar  altitude 
expressed  by: 


bg  = 0.217  - 0.00199  SA  Equation  7-4 

The  fitting  of  this  regression  to  the  skylight  data  is  shown  in  Figure  7-3. 

No  attempt  to  correlate  bg  to  moisture  scale  height  was  made  because  of 
statistical  errors  associated  with  normalization  of  the  data  to  an  equivalent 
no-surround  condition. 

7.3.2  Terrain  Reflectance  Component 

Because  the  vertical  skylight  data  was  altered  specifically  to  eliminate 
surround  reflection  effects,  it  therefore  could  not  be  used  to  formulate  a 
model  of  this  effect.  After  reconstructions  of  the  vertical  spectral  skylight 
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Figure  7-3.  Vertical  Back-Lit  Skylight  Coefficient  Versus  Solar  Altitude 
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samples  were  examined,  it  was  decided  that  this  term  could  satisfactorily 
be  approximated  by  the  one  for  vertical  daylight  irradiance  by  eliminating 
the  equation  change  associate  ' with  specular  reflection  at  low  solar 
altitude  (see  Section  6 of  this  volume). 

7.4  VERY  LOW  SOLAR  ALTITUDE  MODELING 

Solar  altitude  data  as  low  as  8 degrees  was  used  in  developing  the  vertical 
backlit  irradiance  model.  In  addition,  a scaling  coefficient  was  applied 
to  the  8-degree  spectral  curve  to  extend  the  model  from  8 degrees  to  minus 
5 degrees  solar  altitude.  The  scaling  function  was  again  necessary  as  with 
the  other  irradiance  conditions  already  described  because  the  skylight  does 
not  physically  go  to  zero  when  the  equations  do.  The  terrain  reflectance 
term  is  allowed  to  go  to  zero  at  zero  degrees  solar  altitude  as  required 
by  Equation  6-4,  and  it  is  set  to  zero  for  solar  altitudes  less  than  zero 
degrees . 

7.5  MODEL  EQUATIONS  IN  SCAT3 

The  equation  in  final  form  as  modeled  for  use  in  SCAT3  is: 

H^A)  « Ks  • bg(SA)  H$C(A)  cos0  [1-.3  ts(X)m(9)]  tfl(A 

mffll  Equation  7-5 

♦ Ka  • bA(SA)  HSC(X)  cose  tts(X)tfl(X)]  R(A) 


Where:  HySB(X) 


HSC(A) 


ts(X)'ta(X) 


is  the  vertical  spectral  skylight  irradiance, 
in  watts/m2/5  nm. 

is  the  mean  spectral  extraterrestrial  irradiance, 
in  watts/m2/5nm. 

is  the  solar  zenith  angle  (6  also  equals  90° 
minus  the  solar  altitude). 

is  the  vertical  spectral  transmittance  due  to 
scattering  and  absorption. 
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m(e)  is  the  relative  air  mass  corresponding  to  solar 
zenith  angle  0. 

R(X)  is  the  spectral  reflectance  of  the  surrounding 
surface  near  the  receiver. 


^A(SA)  is  the  surround  reflectance  proportionality  constant 
evaluated  for  a specific  solar  altitude,  as  computed 
in  Equation  6-7. 

bg(SA)  is  the  skylight  proportionality  constant  evaluated 
for  a specific  solar  altitude,  as  computed  in 
Equation  7-4,  and 

l 1.0  for  SA  > 6°  or, 

KS  = | p(SA) * for  SA  < 6° 

K = J 1.0,  SA  > 0° 

A \ 0,  SA  < 0°  . 


Typical  vertical  spectral  skylight  j.rradiance  curves  reconstructed  from 
Equation  7-5  are  contained  in  Appendix  C along  with  the  actual  measured 
data. 


* For  a description  of  the  polynomial  p(SA),  see  para.  4.4  and  Appendix  F. 
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APPENDIX  A 

ABSTRACT  DESCRIPTION  OF  ANALYTICAL  SOFTWARE 

A. 1 INTRODUCTION 


The  software  needed  to  analyze  the  radiometric  van  data  was  a combination 
of  package  programs  for  the  IBM  370  and  Hewlett  Packard  9830,  and  software 
written  in-house  for  the  IBM  370,  Hewlett  Packard  9830,  and  Data  General 
Nova. 

In  general,  three  different  types  of  analytical  software  were  used  for  the 
analysis:  (1)  the  transmittance  analysis,  which  used  the  nonlinear  and 
linear  regression  programs  previously  developed  for  the  IBM  370  and  Hewlett 
Packard  9830;  (2)  the  analysis  of  the  scattering  function,  which  used  pro- 
grams developed  for  the  Data  General  Nova  and  a nonlinear  regression  program 
on  the  IBM  370;  and  (3)  the  daylight  horizontal,  skylight  horizontal,  day- 
light vertical,  skylight  vertical, 'and  sky  radiance  analysis,  which  used 
linear  regressions  on  the  IBM  370  and  the  Hewlett  Packard  9830. 

The  equipment  required  in  the  analysis  was:  (1)  an  IBM  370  with  a 3330  disk 
pack,  card  reader,  printer,  an  Eastman  Kodak  Komstar  microfiche  generator, 
and  a 936  Calcomp  plotter;  (2)  a Hewlett  Packard  9830  with  a cassette  unit  and 
plotter;  and  (3)  a Data  General  Nova  with  a nine-track  tape  drive  and  a tele- 
typewriter. 

The  three  different  types  of  analytical  software  involved  are  described  be- 
low, and  a flowchart  that  illustrates  each  one  is  included.  All  of  the  flow- 
charts contain  the  EXTRACT  Program,  which  is  used  to  transfer  data  from  the 
direct  access  files  to  the  sequential  access  files.  (See  Section  7,  Vol  2, 
for  a description  of  the  EXTRACT  Program.) 

A. 2 ANALYTICAL  SOFTWARE  FOR  ANALYSIS  OF  THE  TRANSMITTANCE 

Figure  A-l  is  a flowchart  that  shows  the  programs  used  for  transmittance 
analysis.  Seven  programs  were  involved  in  its  analysis.  The  first  program, 
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EXTRACT,  has  already  been  mentioned. 


The  second  program,  TRALINEAR,  computed  the  coefficients  of  the  linear  re- 
lationship of  airmass  and  log  transmittance  for  each  half  day  in  which  data 
was  collected.  The  coefficients,  slope  and  intercept,  were  then  stored  in 
a disk  file  for  further  processing. 

The  third  program,  CRETSCAT,  has  two  inputs:  (1)  a card  input  that  specifies 
the  1/2  day  for  which  the  analysis  is  to  be  performed,  and  (2)  the  file  that 
was  created  by  TRALINEAR  that  contained  the  coefficients.  The  two  functions 
of  CRETSCAT  are:  (a)  to  eliminate  absorption  and  Rayleigh  scattering  from  the 
slope-estimated  transmittance  so  that  it  has  only  aerosol  scattering  transmit- 
tance (the  slope  term  from  TRALINEAR  was  used  to  estimate  the  transmittance  as 
described  in  Section  2) , and  (b)  to  put  the  data  created  above  in  the  card  format 
which  is  needed  to  run  an  in-house  non-linear  regression  program.  The  output 
from  CRETSCAT  is  a card- formatted  disk  file  of  aerosol  scattering  log  transmit- 
tance, and  a disk  file  containing  the  data  sets  chosen  for  the  analysis. 

The  fourth  program,  TREGANA,  is  a standardized  in-house  nonlinear  regression 
program  that  solves  for  the  coefficients  C and  Ka  referred  to  in  Section  2. 

The  output  from  this  program  consists  of  (1)  the  coefficients  C and  K#  written 
>n  a disk  file  in  card  format,  (2)  a printout  of  the  estimated  value  of  the 
aerosol  scattering  log  transmittance,  and  (3)  a printout  and  plot  of  the  resid- 
uals: that  is,  the  actual  aerosol  scattering  minus  the  estimated. 

The  fifth  program  CRETABSORP  has  two  functions:  (1)  to  isolate  the  absorption 
component  of  the  slope-estimateo  transmittance,  and  (2)  to  create  a card- 
formatted  disk  file  that  is  passed  to  a multiple  linear-regression  program. 

The  absorption  component  is  isolated  by  subtracting  the  Rayleigh  scattering 
and  estimated  aerosol  scattering  from  the  original  slope-estimated  transmittance 
that  was  contained  in  the  file  from  CRETSCAT.  The  disk  file  contains  the  ab- 
sorption component  of  the  transmittance  analysis  plus  data  for  the  water, 
ozone,  and  oxygen  bands.  This  data  is  regressed  agains  the  absorption  com- 
ponent . 


The  sixth  program  is  again  TREGANA.an  in-house  multiple  linear  regression 
in  which  coefficients  are  calculated  for  the  water,  ozone,  and  oxygen  bands 
to  predict  the  absorption  component.  The  input  file  is  from  the  program 
CRETABSORP,  and  the  output  files  consist  of  (1)  the  coefficients  from  the 
regression  which  are  on  a disk  file  that  is  in  card  format,  and  (2)  a print- 
out of  the  estimated  value  of  absorption,  the  residuals  (.the  actual  absorp- 
tion minus  the  estimated),  and  the  statistics  associated  with  the  regression. 

The  last  program  for  the  transmittance  analysis  is  CRETBACKSUB,  which  uses 
(1)  the  coefficients  from  the  non-linear  regression,  (2)  the  coefficients 
from  the  linear  regression,  and  (3)  the  data  and  atmospheric  pressure  that 
were  written  on  the  intercept  and  slope  file  by  CRETSCAT  to  recreate  the  total 
log  transmittance. 

The  final  regressions  between  the  nonlinear  and  linear  coefficients  and 
the  meteorological  parameters  was  done  on  a Hewlett  Packard  9830  calculator. 
These  results  are  discussed  in  Section  2. 

A. 3 ANALYTICAL  SOFTWARE  FOR  THE  SCATTERING  FUNCTION 

The  software  for  this  function  necessitated  use  of  the  Data  General  Nova,  the 
IBM  370,  and  the  Hewlett  Packard  9830.  A flowchart  of  the  Data  General  and 
IBM  370  programs  needed  for  the  analysis  is  shown  in  Figure  A-2. 

The  Data  General  software  was  used  in  creating  the  correction  factors  for 
time  and  path  length.  The  inputs  for  the  Data  General  program  included: 

(1)  a back-up,  nine-track  magnetic  tape  that  contained  irradiance  measurements 
necessary  for  the  time  correction;  and  (2)  the  date,  time,  solar  altitude, 
solar  azimuth,  and  transmittance  that  were  input  manually  at  run  time.  The 
date  and  time  were  used  to  locate  the  data  set  of  interest  on  the  magnetic 
tape;  and  the  solar  altitude,  solar  azimuth,  and  transmittance  were  needed  to 
compute  the  path  length  correction.  When  the  correction  factors  were  calcu- 
lated, they  were  printed  out  on  the  teletype.  After  these  factors  were  checked 
manually  for  validity,  they  were  then  punched  on  cards  for  input  into  an  IBM 
370  program. 


Program  EXTRACT,  which  retrieves  the  required  data  from  the  disk  files  as 
described  above,  places  this  data  on  a sequential  disk  file  ready  to  be 
processed  by  a subsequent  program. 

This  program's  three  main  functions  are  (1)  to  apply  the  correction  factor 
to  the  sky  radiance  data,  (2)  to  extract  data  that  has  an  instrument  elevation 
of  seven  degrees  only,  and  (3)  to  generate  the  disk  file  in  the  card  format 
that  is  necessary  to  operate  the  in-house  non-linear  regression  program. 

This  non-linear  regression  program  determined  the  coefficients  A^,  Aj,  A2>  and 
Aj  of  Equations  3-14  through  3-17  for  each  data  set.  The  output  of  this  pro- 
gram is  a printout  of  the  coefficients  and  a plot  of  the  residuals;  that  is, 
the  actual  data  minus  the  estimated  data. 

A Hewlett  Packard  9830  desk  calculator  was  used  to  complete  the  analysis. 

A linear  regression  program  that  was  written  for  the  calculator  was  used  to 
perform  the  appropriate  regressions. 

A-4 . ANALYTICAL  SOFTWARE  FOR  THE  IRRADIANCE  AND  RADIANCE  DATA 

The  daylight  horizontal,  skylight  horizontal,  daylight  vertical,  skylight 
vertical,  and  sky  radiance  all  use  similar  software.  Figure  A-3  shows  the 
flow  of  the  data  in  the  IBM  computer. 

The  main  purpose  of  the  EXTRACT  Program,  which  was  referred  to  above, is  to 
extract  the  desired  data  from  the  direct  access  files.  Both  irradiometer 
files  were  needed  because  each  irradiometer  was  used  to  record  the  different 
irradiance  types.  If  only  sky  radiance  data  was  processed, only  the  sky  radi- 
ometer file  was  used. 

The  file  was  then  sorted  by  year,  month,  day,  hour,  and  minute  by  an  IBM 
System  sort  routine  and  sent  to  a program  that  converts  the  data  to  even  in- 
crements of  solar  altitude.  To  accomplish  this,  a method  of  linear  inter- 
polation was  used.  The  same  program  also  provides  a time  correction  for  low- 
solar-altitude  data. 
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A permanent  file  is  used  to  store  the  data  that  was  converted  to  even  increments 
of  solar  altitude  and  time- corrected.  This  permanent  file  is  the  data  base 
for  the  analyses  described  in  Sections  3 through  7. 


This  file  was  then  input  to  a program  used  to  prepare  the  data  for  the  regres- 
sion program,  TREGANA.  In  the  case  of  daylight  horizontal  irradiance,  the 
program  subtracted  the  direct  component  from  the  total  measured  irradiance  and 
then  passed  the  data  to  the  linear  regression  program.  In  the  case  of  sky 
radiance,  the  program  converted  the  data  to  log  base  ten  before  passing  data  to 
the  regression  program. 
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APPENDIX  B 

GLOSSARY  OF  VARIABLE  NAMES  IN  SC AT 3 


ABBREVIATIONS 

AM(0) 

AZ 

CATS 

DUST 

EL 

H20(A) 

hd  • [A,  Ta(A),  Ts(A),  SA,  AM(ZEN)] 

Hdh  [A,  SA,  MSH,  R(A),  TA(A) , Tg(A) , 
AM (ZEN)] 


DEFINITION 

Airmass  as  a function  of  solar  zenith 
angle  6. 

Solar  azimuth  angle  in  degrees. 

Angle  between  observer's  ray  and 
sun's  ray. 

Atmospheric  dry  aerosol  content  in 
parts  per  milliliter. 

Terrain  elevation  above  sea  level 
in  kilometers. 

Spectral  structure  for  the  water 
absorption  bands. 

Direct  Solar  Irradiance  Component* 

Daylight  irradiance  on  a horizontal 
surface.* 


Hdv  [A,  SA,  MSH,  R(A),  TA(A),  Tg(A),  Daylight  irradiance  on  a vertical 
^(ZEN)]  surface  perpendicular  to  the  sun's 

azimuth.* 

Hpg  [A,  SA,  MSH,  Ta(A).  Tg(A),  Indirect  downward  scattered  skylight 

AN! (ZEN) ] irradiance  on  a horizontal  surface.* 


Hr  [A,  SA,  MSH,  R(A),  TA(A),  T„(A)  Indirect  downward  scattered  skylight 

AM(ZEN)]  reflected  from  the  terrain. 

Johnson's  extraterrestrial  spectral 
irradiance  outside  the  atmosphere.* 

HgH  [A,  SA,  MSH,  R(A) , TA(A) , Tg(A),  Skylight  irradiance  on  a horizontal 
AM (ZEN)]  surface.* 


HSC(A) 


HgvB  [A,  SA,  MSH,  R(A) , T (A) , Tg(A),  Skylight  irradiance  on  a vertical 
AM(ZEN)]  surface  perpendicular  to  the  sun's 

azimuth  and  pointed  away  from  the 
sun. 

IEL  Instrument  elevation  above  the  horizon 

in  degrees. 

A Wavelength  in  nanometers. 

MSH  Moisture  scale  height  in  centimeters. 

•All  irradiances  have  the  units:  watts/meter2/5  nanometers. 
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APPENDIX  B (Cont'd) 


ABBREVIATIONS 

N[A,  OBS,  CATS,  ZEN,  AM(OBS),  TA(A), 
Tg( A)  MSH,  R(6S0) ] 

o2(a) 

Oj(A) 

OBS 

P 

R(A) 

RH 

SA 

TC 

TF 

Ta  [A,  03(A),  H20(A),  02(X)] 

Ts  [A,  P,  MSH,  DUST] 

ZEN  or  4> 


DEFINITION 

Atmospheric  path  radiance  in  watts/ 
metervsteradian/5  nanometers. 

Spectral  structure  for  the  oxygen 
absorption  band. 

Spectral  structure  for  the  ozone 
absorption  band. 

Angle  between  observer's  line-of- 
sight  and  local  zenith. 

Pressure  in  millibars. 

Spectral  reflectance. 

Relative  humidity  in  percent. 

Solar  altitude  in  degrees. 

Temperature  in  degrees  centigrade. 

Temperature  in  degrees  Fahrenheit 

Spectral  transmittance  due  to  absorp- 
tion. 

Spectral  transmittance  due  to  scattering 
Solar  zenith  angle  in  degrees 


This  program  also  created  the  control  cards  necessary  to  run  the  in-house, 
multiple- linear-regression  program  TREGANA.  This  program  then:  (1)  com- 
puted the  appropriate  coefficients,  (2)  printed  and  card  punched  the  co- 
efficients, and  (3)  printed  the  statistics  of  the  regression. 

A subsequent  program  then  performed  the  calculations  needed  to  run  a second 
regression,  which  included  calculating  the  residuals  between  the  results  of 
the  first  regression  and  the  data  that  went  into  the  first  regression.  For 
the  irradiance  data,  only  the  days  with  snow  albedo  were  passed  to  the  second 
regression;  this  program  therefore  extracted  only  those  days  with  this  type 
of  albedo.  The  coefficients  from  the  linear  regressions  where  then  regressed 
against  solar  altitude  and/or  moisture  scale  height.  For  this  purpose,  a 
Hewlett  Packard  9830  calculator  multiple  linear  regression  program  was  used. 


Several  computers  were  used  for  different  purposes;  for  example,  the  IBM  370 
was  used  for  large-scale  data  analysis,  the  Hewlett  Packard  9830  was  used  for 
small-scale  data  analysis  and  testing,  and  the  Data  General  was  used  to  re- 
format data  to  make  them  compatible  with  the  IBM  370  system.  Because  of  the 
large  volume  of  printout  from  the  IBM  most  of  it  was  written  to  microfiche 
to  provide  easy  accessibility. 


APPENDIX  C 


EXAMPLES  OF  SCAT3  RECONSTRUCTION  OF  ACTUAL  DATA 


Estimates  of  actual  measured  data  for  a selected  cross  section  of 
meteorological  and  solar  altitude  conditions  are  shown  in  Figures  C-l 
through  C-84.  Reconstruction  was  accomplished  using  the  surface  meteo- 
rology algorithms  of  SCAT3,  substituting  in  the  exact  measured  pressure 
temperature,  relative  humidity,  and  geometry  given  in  Table  C-l  below. 
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INPUT  PARAMETER  SUfWARY  FOR  SCAT3  RECONSTRUCTION  EXAMPLES 
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APPENDIX  D 


SCAT3  MODEL  RECONSTRUCTED  SPECTRAL  TRANSMITTANCE 


The  following  Figures  D-l  through  D-3  illustrate  the  dependency  of  the 
SCAT3  model  for  spectral  transmittance  on  surface  pressure,  relative  humidity, 
and  temperature. 
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APPENDIX  E 

TIME-FOR-MEASUREMENT  CORRECTION  REQUIRED 
FOR  LOW  SOLAR  ALTITUDE  SKY  RADIANCE  DATA 


Because  sky  conditions  change  rapidly  at  low  solar  altitudes,  the 
measurements  made  at  these  altitudes  need  to  be  corrected  by  use  of 
a time-correction  factor.  For  this  purpose,  a correction  factor  was 
determined  from  measurements  taken  at  3S0  nm  immediately  following  a 
full  spectral  sequence  of  measurements  that  compared  these  readings 
with  the  initial  spectral  measurement,  which  was  at  350  nm. 

The  percentage  change  was  divided  equally  into  173  parts  since  the 
first  measurement  was  taken  at  a starting  time  of  1 unit  and  the  last 
measurement  was  taken  at  a starting  time  of  173  units. 

The  final  equation  for  correcting  data  was: 

DataCorr(>)  ’ Data  Uticorr(i)  [1  * 1 


Where:  DataCorr(X) 


is  the  time- corrected  spectral  measure- 
ment 


DataUncorr 

n 


(X)  is  the  uncorrected  spectral  measure- 
ment 

is  the  nth  spectral  sample 


173  is  the  number  of  spectral  samples  in- 

cluding the  first  350  nm  sample  and  the 
last  350  nm  spectral  sample,  and 

pr  _ Data(l)  - Data Cl 73) 

Data(l) 


Not  all  spectral  sequences  at  solar  altitudes  below  8°  had  a 350  nm  mea- 
surement taken  immediatel y following.  In  order  to  correct  those  data 
that  did  not  have  the  extra  measurement,  a function  for  PC  was  formulated 
from  the  sequences  that  had  had  the  extra  measurement.  PC  was  found  to 
be  a function  of  solar  altitude. 


E-l 


Thus,  from  repeated  measurements  at  various  solar  altitudes  for  the 
daylight  vertical,  skylight  vertical,  daylight  horizontal,  and  skylight 
horizontal  irradiances  the  following  equation  was  determined  for  the 
percentage  change  as  a function  of  solar  altitude: 

PC  - .23784  - .041119  • SA  ♦ .0025274  • SA2 


where  SA  is  between  -1°  and  8®. 


APPENDIX  F 


VERY  LOW  SOLAR  ALTITUDE  SCALING  OF  SPECTRAL  QUANTITIES 

At  solar  altitudes  of  eight  degrees  or  below,  either  sufficient  data 
was  not  always  available  to  test  the  model  or  the  model  was  not  ade- 
quately describing  the  data,  whether  it  was  radiance  or  irradiance  da- 
ta. At  lower  solar  altitudes  than  those  thought  to  be  the  minimum  so- 
lar altitude  (MSA)  for  an  accurate  fit,  the  model  was  extrapolated  down 
to  -5  degrees  by  multiplying  the  model  results  at  the  MSA  by  a scaling 
coefficient.  This  scaling  coefficient  is  solar- altitude -dependent  and 
is  based  on  horizontal  daylight  irradiance.  The  scaling  coefficient 
(SC)  is  considered  valid  down  to  -5  degrees  solar  altitude.  This  co- 
efficient is: 

f (SA) 
in  1 

SC  ■ — — Equation  F-l 

10fl<MSA) 

The  factors  fj(SA)  and  fj(MSA)  are  computed  using  the  following  polynomial: 

f I (SA  or  MSA)  ■ 2.112  + .1575  SA  - .00721  SA2  Equation  F-2 

where  SA  is  the  solar  altitude  between  -5*  and  90°,  and  MSA  is  the 
lowest,  valid  solar  altitude  between  -5®  and  90°. 


Application  of  the  coefficient  to  haze  radiance  or  irradiance  gives  the 
following  equation: 


H(X,SA)  - SC  • H(A,  MSA)  Equation  F-3 


F-l 


Each  irradiance  type  and  the  haze  radiance  have  a different  minimum 
valid  solar  altitude  (MSA).  For  this  reason,  when  the  MSA  is  involved, 
reference  should  be  made  to  the  different  sections  of  this  report  to 
find  the  MSA  associated  with  each  irradiance  or  radiance. 


F-2 


APPENDIX  G 

CALCULATION  OF  APPARENT  SOLAR  ALTITUDE 


The  continuous  change  in  the  index  of  refraction  of  the  atmosphere  as  a 
function  of  altitude  causes  the  apparent  position  of  the  sun  to  be  higher 
than  its  actual  geocentric  position.  For  a typical  atmospheric  density 
profile,  this  solar  altitude  error  was  computed  by  McClatchey  et  al*  for 


solar  zenith  angles  in  excess  of 

75°.  These  data  are  reproduced  below. 

Apparent 

Solar  Zenith  Angle 

Geocentric 

Solar  Zenith  Angle 

(in  Degrees) 

(In  Degrees) 

75 

75.061 

80 

80.090 

83 

83.124 

85 

85.167 

86 

86.199 

87 

87.248 

88 

88.311 

88.5 

88.858 

89 

89.417 

89.5 

89.997 

90 

90.570 

Since  the  SCAT3  atmospheric  model 

uses  trigonometric  functions  to  compute 

the  parametric  spectral  vectors  for  several  radiometric  quantities,  it  was 

necessary  to  substitute  apparent 

solar  altitude  in  the  equations  rather  than 

true,  geocentric  solar  altitude. 

The  conversion  between  geocentric  and 

apparent  solar  altitude  is  done  in  the  program  by  using  the  following 

polynomial  equation  : 

SA  - SA  ♦ .544  - 

a g 

.1572.SA„  ♦ .01764*SA2b 
g * 

* McClatchey  et  al,  "Optical  Properties  of  the  Atmosphere,"  Air  Force 
;f  Cambridge  Research  Labs  71-0279,  May  1971,  Table  6,  p 41. 

G-l 


Where: 


is  the  geocentric  solar  altitude  in  degrees,  and 


' SA  is  the  apparent  solar  altitude  in  degrees. 

Si 

This  equation  was  derived  from  the  above  table  data  by  first  converting 
the  data  from  zenith  to  altitude  angles  (90°  - solar  zenith  angle)  and  then 
regressing  the  apparent  solar  altitudes,  SA  , against  the  geocentric  solar 

o 

altitudes,  SA^. 

No  correction  is  applied  to  solar  altitudes  higher  than  5°.  Below  0°  apparent 
solar  altitude,  corrected  values  are  computed,  but  they  have  no  significance 
for  the  SCAT3  equations. 
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APPENDIX  H 


1 


PATH-LENGTH  CORRECTION  FOR  THE  POLAR  SCATTER  FUNCTION 

Each  particular  volume  of  air  in  the  sky  that  is  observed  by  the  sky 
radiometer  receives  a different  amount  of  light  depending  on  the  extinc- 
tion coefficient  of  the  atmosphere,  and  the  atmospheric  distance  each 
point  is  away  from  the  sun.  In  addition,  the  sky  spectroradiometer  re- 
ceives an  amount  of  light  from  this  volume  of  air  along  its  line-of- 
sight  that  is  diminished  by  the  air  path  length  to  the  radiometer.  As  a 
result,  the  total  relative  amount  of  light  reaching  the  radiometer  can  be 
expressed  in  terms  of  the  following  quantities: 

a.  The  amount  of  light  reaching  a particular  volume  of  air  whose 
position  lies  along  the  observer's  line-of-sight, 

b.  The  amount  of  light  reaching  the  observer  from  that  volume  of 
air,  and 

c.  The  scattering  of  light  by  the  volume  of  air  toward  the  obser- 
ver. 

The  geometry  involved  in  the  above  is  shown  in  Figure  H-l  for  the  volume 
of  air  at  point  P in  the  atmosphere.  Two  integrations  are  involved:  (1) 
over  the  path  defined  by  the  sun-to-volume  line,  and  (2)  the  volume-to- 
radiometer  path  defined  by  the  observer's  line-of-sight.  The  integral  that 
expresses  energy  reaching  the  observer  (or  radiometer)  is: 

Lmax 

PL(AZ)  « J Tp_0*  Tg_p(AZ)  • Hgc  • S • dL  Equation  H-l 

Where:  PL(AZ)  is  the  relative  amount  of  light  reaching  the 
observer  at  point  0 from  a volume  of  air  at 
point  P and  for  an  observer's  azimuth  AZ  rela- 
tive to  the  sun, 

T£  p(AZ)  is  the  transmittance  of  light  between  the  edge 
of  the  atmosphere,  point  E,  to  a volume  of  air 
along  the  observer's  line-of-sight,  point  P, 
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Figure  H-l.  Definition  of  Geometry  Conditions  of  Integration 
of  the  Sampled  Volume  of  Atmosphere 


w* 


is  the  transmittance  of  light  between  the 
volume  of  air  at  point  P along  the  observer's 
line-of-sight  to  the  observer  at  point  0, 

is  the  distance  along  the  line-of-sight  of  the 
observer  (or  instrument  field  of  view)  with 
zero  corresponding  to  the  earth's  surface  and 
Lmax  to  the  intersection  of  the  line  with  the 
effective  edge  of  the  atmosphere, 

is  the  solar  irradiance  outside  the  atmosphere, 
and 

is  the  fraction  of  light  scattered  by  the  volume 
of  air  along  the  observer's  line-of-sight  toward 
the  observer.  (S  is  taken  to  be  proportional 
to  the  atmospheric  density  of  the  volume  of  air 
at  height  h above  the  earth  and  can  be  expressed 
as  S = e~ ■ 1535h) . 


The  following  equation  describes  the  transmittance  of  light  reaching  a vol- 
ume of  air  located  along  the  observer's  line-of-sight  as  a function  of  alti- 
tude and  path  length.  In  Figure  H-2,  the  distance  r corresponds  to  this 


transmittance. 


-rmax 

J ane' 

*'r\  O 


1535H • 


Equation  H-2 


Where: 


is  the  distance  along  a line  between  the  sun  and  the 
volume  of  atmosphere  under  consideration  as  measured 
from  the  edge  of  the  atmosphere. 


a is  the  extinction  coefficient  of  the  atmosphere  at  the 

earth's  surface  calculated  for  each  day,* * 

rmax  is  the  distance  from  the  outer  edge  of  the  atmosphere 
to  the  volume  of  atmosphere  under  consideration,  and 

H is  the  perpendicular  height  above  the  earth  to  point  r. 


* The  daily  calculation  of  atmospheric  transmittance  is  explained  in 

Section  2 above. 
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Figure  H-2.  Trigonometric  Definitions  for  an  Atmospheric 
Volume  at  Point  P 


The  equation  for  rmax,  the  distance  between  point  E and  point  P,  is 
derived  from  the  Pythogorean  theorem  (see  Figure  H-2)  and  is  expressed 
algebraically  as  follows: 


rmax  * |(6400>h)  • cosZ  i / 64302  - (6400«-h)2  sinzZ  Equation  H-3 


Where:  h is  the  height  of  the  volume  of  air  above  the 

earth's  surface, 

Z is  the  solar  zenith  angle  at  point  L along  the 
observer's  line-of-sight. 


6400  is  the  radius  of  the  earth  in  kilometers,  and 


30  is  the  effective  height  of  the  atmosphere  in 
kilometers. 


The  equation  for  H is  derived  using  the  law  of  cosines  (see  Figure  H-2)  and 
becomes : 

h ■ / (6400+h)2  ♦ (rmax-r)2  ♦ 2(6400+h) (rmax-r)  cosZ  - 6400  (Equation  H-4) 

The  height  h is  also  derived  using  the  law  of  cosines  (see  Figure  H-3  ) and 
is  expressed  algebraically  as  follows: 


h - ✓ 64002  ♦ s2  - cos(ALS)  • s • 6400  - 6400  Equation  H-5 


Where:  s is  the  distance  from  the  observer  to  a point  along 
the  line-of-sight,  and 

ALS  is  90°  plus  the  angle  between  the  observer's 
line-of-sight  and  the  local  horizon. 


The  cosine  of  the  zenith  angle,  Z,  is  derived  from  the  spherical  law  of 
cosines  (see  Figure  H-4);  it  is  expressed  algebraically  as  follows: 

CosZ  ■ cosZ  • cos$  ♦ sinZQ  • sin$  • cos0  Equation  H-6 


H-5 


Figure  H-3.  Trigonometric  Definition  of  the  Angle  $ 
and  Altitude  h 


Figure  H-4 . Spherical  Trigonometric  Definition  of  Angle  Z 


Where:  2 is  the  observer's  zenith  angle, 

♦ is  the  arc  between  the  observer  and  the  nadir  of 
a point  along  the  observer's  line-of-sight,  and 

6 is  the  solar  azimuth  relative  to  the  azimuth  of  the 
observer’s  line-of-sight  (90°  - AZ). 

The  equation  for  sin$  is  derived  using  the  law  of  sines  (see  Figure  H-3) 


s • sin(ALS) 
6400+h 


Equation  H-7 


After  the  relative  amount  of  light  reaching  each  volume  along  the  observer's 
lxne-of-sight  is  determined,  the  amount  of  light  reaching  the  observer  after 
it  is  scattered  by  each  volume  must  be  considered.  The  transmittance  of  light 
between  the  volume  of  air  and  the  observer  is  defined  in  a manner  similar  to 
that  resulting  in  Equation  H-2,  and  becomes: 


\ - . 1535h , 

V ds 


Equation  H-8 


Where : 


smax  is  the  distance  between  the  observer  and  the  edge 
of  atmosphere  in  the  direction  of  the  observer's 
line-of-sight,  and 

s is  the  increment  along  the  line  from  the  volume  of 
air  to  the  observer. 


The  relative  amount  of  light  reaching  the  observer  at  the  observer's  azimuth 
relative  to  the  sun,  AZ,  is  determined  from  Equation  H-l.  it  is  used  to 
normalize  the  sky  radiance  to  one  point  in  the  sky.  This  normalizing  is  done 
by  evaluating  PL(AZ)  in  Equation  H-l  at  the  azimuth  corresponding  to  direct 
backscatter  (away  from  the  sun)  and  dividing  through  by  PL(AZ)  evaluated  at 
each  of  the  other  azimuths  of  measurement.  As  indicated  below  and  in  para- 
graph 3.4  above,  the  corrected  path  radiance  is: 


N" (AZ)  • N ' (AZ)  • 


m 


Where: 


N"(AZ)  is  the  path  length  and  time- corrected  sky  radiance 
at  azimuth  angle  AZ, 

N'(AZ)  is  the  sky  radiance  at  azimuth  angle  AZ  before 
path- length  correction,  and 

AZ  is  the  instrument's  azimuth  angle  relative  to  the 
solar  azimuth.  (The  direction  away  from  the  sun 
is  0‘). 
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APPENDIX  I 

CORRECTION  OF  MEAN  EXTRA-TERRESTRIAL  IRRADIANCE  BY  TIME-OF-YEAR 

The  SCAT3  model  incorporates  a function  to  change  the  extra  terrestrial 
spectral  irradiance  with  time-of-year  according  to  the  orbital  eccentricity 
of  the  earth.  This  correction  is  described  below. 

The  values  for  the  Solar  Constant  are  the  mean  of  the  extraterrestrial  ir- 
radiance values,  and  are  given  for  the  mean  distance  between  the  sun  and 
earth  (one  astronomical  unit  distance) . The  annual  variation  in  irradiance 
due  to  the  cyclical  variation  in  earth-sun  distance  amounts  to  approximately 
7%,  and  can  be  determined  with  much  greater  accuracy  than  the  Solar  Constant 
itself.*  This  variation  can  be  computed  directly  from  an  eccentricity  of 
the  earth's  orbit  and  from  the  mean  anomaly  (the  angular  difference  between 
the  mean  solar  longitude  and  the  mean  solar  longitude  of  perigee) . The  com- 
putational procedure  as  used  in  SCAT3  is  outlined  below. 

The  irradiance  factor  applied  to  the  mean  iTradiance  is  given  by: 

K * (I/R)2  Equation  I- 

where  R is  the  earth-sun  radius  vector.  This  radius  vector  (R)  is  computed 
from: 

R = 1 ♦ (1/2)  e2  - (e  cos  M)  -[(1/2)  e2  cos2  M] 
where  e is  the  eccentricity  of  the  earth's  orbit,  and  M is  the  mean  anomoly. 

The  eccentricity,  which  diminishes  constantly,  can  be  computed  from: 

e * .01675  - .00004  T 

where  T is  the  number  of  centuries  elapsed  from  the  epoch  beginning  noon, 

1 January  1900,  each  century  being  of  36525  ephemeris  days  (Julian  Calendar). 

Therefore,  T » d/36525. 

* Explanatory  Supplement  to  the  Astronomical  Ephemeris  and  the  American 

Ephemeris  and  Nautical  Almanac,  H.M.  Nautical  Almanac  Office,  Her  Majesty's 
Stationary  Office,  1961. 


The  mean  anomaly  is  computed  from: 


M « 358.47584S  ♦ 0.9858600267d  - 0.0000112  D2  - 0.00000007  D3 

Where:  d is  the  number  of  ephemeris  days  elapsed  from  noon,  1 January 
1900,  and 
D is  d/10000. 


